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POUKWOltD 


ThlR report documontfl tho work conductor! under Coni mot NAH8-2708G, 
"Htudy of Cnvltntlng Inducer lnfltnbll!tlon," during tho period 1 .July 1071 through 
11} May 1078. Tho work was flponnorod by tho George (’. Marshall Spoor' Flight 
Center, National Aoronautlcn and Hjiaoe Administration, Mnrflhnll Space Flight 
Center, Alabama, and wan nrlnilnlnterod technically by Mr, II, P. Stlnnon, .Tr. 

Pratt & Whitney Aircraft '« Florida UoHoaroh and Development Center at 
Wont Pnlm Dench, Florida, wan tho contractor, and Mr. W. K, Young wan tho 
Program Manager. All work wan performed at Fit DC with the exception of the 
cuBcudo tests, which wore conducted at United Aircraft Research Laboratories, 
Fast Hurtford, Connecticut, by Mr, W. F. Taylor. Inducer tests woro conductorl 
by Mr, A. F. Wemmell. 


ADSTRACT 


An analytic and experimental Investigation into the causes and mechanisms 
of cavltatlng inducer instabilities was conducted. Hydrofoil cascado tests wore 
performed, during which cavity sizes were measured. The measured data wer > 
used, along with Inducer data and potential flow predictions, to refine an analysis 
for the prediction of Inducer blade suction surface cavitation cavity volume. 

Cavity volume predictions were incorporated Into a linearized system model, 
and Instability predictions for an Inducer water test loop were generated. Inducer 
tests were conducted and Instability predictions correlated favorably with measured 
Instability data. 


SUMMARY 


An anatytlc and experimental Investigation Into the causes and mechunlmns 
of cavltatlng Inducer Instabilities was conducted. Two possible Instability mech- 
anisms were analytically identified: (1) a "continuity" mechanism that results 
from the response of Inducer cavitation volume to Inlet flowrate; and (2) a "per- 
formance" mechanism that results from the relationship of inducer headrlse to 
Inlet flowrate In the cavltatlng head falloff region of operation. The continuity 
mechanism requires only the presence of flowrate-sensltlvo Inducer cavltat'on 
for an Instability to be possible, and this mechanism can cause tnstubltltles at 
relatively high cavitation numbers where there Is no head falloff (l.e. , In the 
usual rango of Inducer operation). Tho performance mechnnlBtn Is operable 
In the hoad falloff region. An analytic system model, In which the Inducer and 
Its system are treated us u series of lumped resistances, Inortances, and 
compliances, can predict tho occummco of both types of Instabilities. The 
model consists of a series of equations that define (1) the Inlet and dtschai’go 
line: resistances, Inertnnces, and compliances; and (2) the Inducer: cavltatlng 
hoad vs Inlet pressure and flow curves (heud falloff map), head vs flow curves 
(head flow map), and cavitation volume vs pressure and flow curves (cavitation 
votumo map). The rati of change of cavitation volume with Inlet pressure Is 
dcflnod ns "pressure compliance" and the rate of change of cnvltutlon volume 


with flowrate Ir do f I nod hr "flow compliance, " Tho nioflt algnlfloant pnramotopR, 
nnd tho moot difficult to accurately doflno with regard to "continuity" insta- 
bllltlofl, nro proRRuro oompllnnco nnd flow campllnnoo, "Performance" typo 
InntnhllltloR require definition of tho bond fulloff map, which cun uminllv be 
obtained through tho uoo of para metric moanurnd data from almllar Inducers. 

Two oxporlmontn wore conducted: (J) a cancado experiment whono 
objective wao to obtain moanurnd cavity geometry data for refinement o' Induce r 
blade Ruction surface cavitation volume prediction (pros nu re compliance and 
flow compliance)} and (8) an Inducer experiment whose objective was to provide 
data for substantiation of the analytic Instability predictions technique. Results 
of the cascade experiment wore largely Inconclusive because of apparent flew 
nonunlformltlcs. Cavity volume predictions, therefore, wore correlated and 
refined using limited Inducer measured data and potential flow predictions. 

Tho Inducor oxporlmont resulted In tho measurement of "continuity" and "per- 
formance" Instabilities whoso characteristics correlated with llnoar Inducor 
system model predictions In all significant areas. Agree mont between pre- 
dicted nnd measured frequoney magnitude was poor, but measured frequencies 
were known to bo affected by system compliances that were not accounted for 
In tho prediction unulysls. Such compliances were of sufficient magnitude to 
explain the difference In frequencies. The analytic model of the Inducor system 
Indicated that "flow compliance" was the most significant parameter with regard 
to system stability. If tho Inducer had no flowrato sensitive cavitation, It would 
have boon predicted to bo stable regardless of the value of other parameters. 

It was concluded that the analytically derived Instability mechanisms 
explain the Instabilities encountered In the Inducer test program. The analytic 
technique employed can be generally used to predict Instabilities and define 
stabilizing system changes. 


SECTION 1 
INTRODUCTION 


Inducers lire widely lined In liquid rocket engine propellant feed n.vnleniH to 
fnereane propellant prennure before the propellant entern the main turliopumpn, 
The turho|Himpn are thereby permitted to operate at relatively high rotational 
opcode with nufflelent Inlet prennure to preclude cavitation performance loan, 

The Indueern themnelven operate with nlgnlfleanl cavitation, but they are denlgned 
Much that the effect of cavllatlon en overall InduceiHurbopump nyntem prennure 
rise and efficiency are minimal. 

The exlulence of cavitation In Die feed Hyntem ralncn the ponnlblllty that 
"Hell- Induced InutubllltloH" may occur. Thone Instabilities ure churn cuu'l/ed by 
oHctllatlonn In Inducer Inlet and dlneh irge prennure which can, If their amplitude 
Ih Hufflclently large, rcnult In unaocoptable varlullonH In engine t brunt becauHe 
of tholr effect on propellant llowrate. Such oHdllatloiiH have been ubnerved 
during Inducer and vehicle Hyntem Ic4h by a number of InveHtlgatorn, and their 
occurrence has been generally linked to Inducer blade Huctlon surface cavitation; 
however, there has been little huocobs In developing an analytic method for 
predicting the occurrence of on tmitablllty or for the definition of system changes 
that would be required to avoid instabilities. 

The objective of this program was to develop an analytical system capable 
of predicting the occurrence ot sell-induced Instabilities In cavltatlng Inducers. 
The technical effort consisted of four phases, the scope of which can he described 
as follows: 

t, Review of Existing Models - Literature concerning Instability 

models was reviewed and the most promising approach selected, 

II, Design of Experiments - Two experiments were designed: A 

two-dimensional cascade experiment to provide measured cavi- 
tation cavity geometry data as a function of Incidence angle and 
Inlet pressure, and rotating Inducer experiment to provide 
measured instability data over a range of operating points, 

ill. Conduct Experiments - The experiments planned In Phase II 
were conducted. 

IV. Data Correlation - Cavity geometry data from the cascade 

tests was correlated with analytic predictions from the cavity 
model of (1)*, and the model was refined as Indicated. The 
developed cavity model was used, in conjunction with a 
system model, to predict the unstable operating region of an 
Inducer. 

The program results Indicate that nn analytic system modeling technique 
that accounts for (1) Inducer Inlet and discharge lino characteristics, (2) inducer 
head rise as a function of cavitation number and flow coefficient, and (3) Inducer 


♦Underlined numbers In parentheses denote references, which arc listed In 
Section 7 on page 143. 





3 


blade auction nurfaeo cavity volume no a function of cavitation number and flow 
coefficient can predict the occurrence of an (notability* and the technique can bo 
unnd to define atabUUlng ayatoro chiuipon, Tleaultn are reported In detail In 
the following noetlonn. 


SECTION 2 

REVIEW OF PREVIOUS WORK 


A lltornturo survey was conducted to ldontlfy reported Instances of 
Instability and related modeling attomptn. Literature pertinent to both rotating 
inducer and stationary tunnel testa was reviewed because apparently similar 
Instabilities have boon noted In both. Because of flow similarities between 
inducor and cascade flows, It aoomod llkoly that a modeling approach that would 
oxplnln Inducor Instabilities would also bo adaptable to two-dimensional water 
tunnel Instabilities. Verification of this similarity would demonstrate that the 
fluid dynamic offoets that load to an Instability nro understood and can bo 
adequately modeled, Increasing the degree of confidence In the Inducor Instability 
model . 

The literature that was found In the survey Is listed In References (2) 
through (22), The reported experimental observations and modeling attempts 
are summarized In the following sections, after which our conclusions arc 
discussed. 

2. 1 EXPERIMENTAL OBSERVATIONS 

Cavltatlon-induced oscillations in water tunnel tests of a single hydrofoil 
were reported by Wade and Acosta (2) who studied the growth of cavitation on a 
plano-convex hydrofoil as Inlet pressure was lowered. They were the first In- 
vestigators to report frequencies and amplitudes during unstable cavitation on 
a hydrofoil. They observed that the cavitation cavity, which began at the hydro- 
foil leading edge on the suction side, was steady and that measured hydrofoil 
normal force amplitudes were small whenever cavity length was less than about 
60% of chord length. With further lowering of Inlet pressure the characteristic 
low frequency, large amplitude cavltatlon-induced oscillations began. With still 
further lowering of Inlet pressure the oscillations persisted, with hydrofoil force 
and cavity length amplitudes first increasing and then decreasing, until the cavity 
reached a length of about 120% chord. At lower Inlet pressures and longer cavities 
the flow became quite steady. Cavity length oscillations reached a peak-to-peak 
amplitude of about 1/2 chord, with hydrofoil force varying ±10% of mean value 
in phase with the cavity. The Instabilities occurred In a region where hydrofoil 
mean lift and drag forces first Increased, then decreased as Inlet pressure was 
lowered. Typical frequencies found were about 12 to 26 Hz, depending upon 
tunnel speed, angle of attack, and cavitation number. Two regions of high fre- 
quency (270 Hz and 50-60 Hz), low amplitude "noise" were also found, but these 
do not appear to have a direct bearing on the Instability presently In question. 
During one cycle of oscillation, the cavity was seen to grow smoothly from Its 
minimum length until, as It approached the hydrofoil trailing edge, a reentrant 
Jet formed and began to gradually fill the rearward portion of the cavity. The 
cavity surface became Irregular and a large volume of cavity was abruptly shed 
Into the stream. The cycle then repeated. No mention is made of tip vortex 
cavitation, nor Is any appreciable amount discernible in the photographs. 

Kaplan and Lehman (3) Investigated cavitation on a semiwedge hydrofoil 
mounted In a 2-degree of freedom suspension In a water tunnel, finding that 
"for cavity lengths In the range from 0.5 to 1. 0 chord the cavity Itself becomes 
unstable and severe fluctuations occur In the entire system. " The surface 
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cnvity began at tho hydrofoil lending odgo. Tho hydrofoil wna cantilevered from 
ono end (essentially flush with tunnel wall) while the other end protruded Into tho 
froostroam. Tip vortex envltntton oocurrod on tho free end, nnd "—apponrfl to bo 
lndopondont uf the envlty on tho foil upper surface. » A trneo of hydrofoil rotn- 
tlonnl and translational motion at ono operating point (cavity length of 70% chord) 
Indicator* a frequency of oscillation of 8 II/. No data are presented that would 
allow a comparison of tho unstable region with tho hydrofoil lift breakdown point. 

Bosch (4) tostod a somlwodgo hydrofoil mounted In a 2-dogroo of froodom 
suspension In a water tunnol, finding that oscillations In oavlty length, with nmpll- 
udes ns much as 25% chord, occurred when tho cnvity was botweon 40 nnd 140% 
chord length. Frequency of oscillation varlod from 4 to 18 Hz, docroaslng ns 
cavity longth Increased. Tho cavity sprang from tho hydrofoil loading odgo. 

No mention Is mado of Up vortex cavitation, nor of how tho unstable region com- 
pared to the hydrofoil lift bronkdown point. H 


Wade and Acosta (5) experimented with two-dimensional cascades of 3 to 
5 plnno-convox hydrofoils In a water tunnel and found an instability similar to 
that reported for the single hydrofoil of the same shapo (2). They were the first 
to experiment with a cascade far enough Into the cavltatlng region to uncover 
unstable cavitation. When the cavities, which extended from the blade leading 
edges, became about 1/3 chord In length, periodic oscillations of the cavities 
and static pressure upstream of the cascade occurred. With further reduction 
of tunnel pressure the cavities became longer than the chord, and oscillations 
ceased. Static pressure oscillations reached an amplitude of about 1/2 the 
approaching velocity head at a frequency of about 12 Hz. Wade and Acosta note 
that there was " — some time lag in the development of oscillating cavitation 
from one vane to the next — . " Presumably this means that the cavities on some 
hydrofoils began to oscillate at higher values of tunnel pressure than the others. 

It is not stated whether the cavities oscillated In phase once they all began to 
oscillate. The only reference to tip vortex cavitation Is a statement that 
— noticeable tip clearance cavitation — " occurred. No mention Is made of a 
propagating cavitation pattern. The data Indicate that the Instabilities occurred 
at or after cascade lift breakdown began. 


Wade and Acosta (6) Investigated another cavltatlng cascade to study 
Instabilities. They measured the amplitudes and phase angles of the cascade 
lift and drag forces and of pressures and velocities upstream and downstream 
ot the cascade. Motion pictures were made of the cavity oscillations. Both a 
surface cavity springing from the leading edge and tip vortex cavities at each 
end of the hydrofoil were formed. Oscillations began when the cavities became 
approximately 60% chord In length, and affected all measured parameters. "This 
oscillating or chugging mode of cavitation persists and grows more Intense as the 
pressure Is lowered until the maximum length of the cavity bubble formed during 
the oscillation cycle Is approximately 16 to 20% longer than the chord of the g 
hydrofoil. At this time the fluctuations In lift and pressure throughout the 
cascade are greatest. Then with a slight, but still further decrease In pressure, 
the cavities become considerably longer than the chord and the oscillations cease 
entirely—. " Frequencies found were from 8 to 20 Hz. From high speed motion 
pictures an attempt W as made to detect any propagating phenomenon that might 
be responsible for the observed Instabilities, The authors state after many 
repeated viewings of these motion pictures we could not conclude that there was 
a progressive disturbance across the cascade In the sense of a propagating stall 
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which Ih responsible for the Instability thnt Is obsorvod. " Tho datn Indlcato that 
tho Instabilities occurred In a band around tho point whoro cnscado lift breakdown 
occur b. 

Taylor, Murrln and Columbo (7) Investigated tho porformanco and point 
of Incipient cavitation of a cascade of double circular-arc hydrofoils, and notod 
a propagating form of cavitation, Photographs show a slnglo Irregular-shaped 
surfaco cavity propagating along tho cascndo. Propagation voloelty was approxi- 
mately .') m/s (10 ft/soc). Thoro was no visible tip vortox cavitation. The 
authors attribute the phenomenon to nonuniform Inlot conditions, and states "when 
the Inlet flow nonuniformity was loss sovox’o — , the cavitation did not propagato 
and tho bubble size changed simultaneously on all hydrofoils.'* Tho roforonco 
to simultaneous chango rofers to change when Inlet prossuro was changed, slnco 
tho Investigation did not go sufficiently far Into the cavltntlng region for a severe 
cavltatlon-lnduced oscillation to havo occurred. 

Acosta (8) reported cavity oscillations on a four-bladed helical Inducer 
with radial loading edges. Ho found that as Inlet pressure was lowered a "patch 
of cavitation" with a "frosty appearance, " formed at the blade tip. Photographs 
show tho cavitation originating near the leading edge on the blade suction side. 
Acosta states that the "—greatest part of the fuzzy cavitation patch arises from 
a tip clearance flow — " and " — Is confined largely to the outer portions of the 
annulus, but at the lowest cavitation numbers It does occur from root to tip. " 

At low Inlet pressure the cavities formed only on alternate blades, but not always 
on the same blades. At high flowrates the alternate blade cavitation pattern was 
stable. At flowrates lower than the best efficiency point the alternate blade 
cavitation appeared to propagate from blade to blade, with the frequency of 
propagation decreasing as Inlet pressure was lowered and decreasing to zero 
just before head breakdown of the Inducer. No numerical oscillation amplitudes 
or frequencies were reported. Three modifications to the Inducer were made in 
an attempt to suppress the oscillations. Increasing the tip clearance offered 
some help; sweeping back the blade leading edge depressed the occurrence to 
lower cavitation numbers; and using a variable head helix (0. 1 to 0. 16 rad, 

6 to 9 deg) greatly reduced the extent and severity of the oscillating mode. In a 
discussion of (8) by Iura, observations are presented of alternate blade cavities 
on four-bladed Inducers which began to oscillate in length as cavitation number 
was lowered. In one Instance a rotating propagation pattern at one-tenth rotor 
speed was observed. The oscillations are reported to have ceased once head 
breakdown occurred. 

Wood (9) visually noted cavity oscillations In three different mixed flow 
Impellers. In one, a three-bladed model, a low frequency (18 Hz) oscillation 
of tho cavity formed at the Impeller inlet occurred. In the other two Impellers, 
with four blades and eight blades each, higher frequency 142 to 176 Hz) oscilla- 
tions occurred In cavitation formations In the rear channels. The rear channel 
cavitation extended across much of the channel, for some conditions completely 
filling the channel, and appeared to be unrelated to the "leading edge cavitation 
formations. " Wood states that rear channel cavitation was "Instigated by the 
tip cavitation vortexes. " No description of the loading edge cavitation is given. 
Wood docs not mention alternate blade cavitation nor circumferential propagation, 
but notes that he had a limited field of vision. It may be of significance that the 
low frequency Impeller Inlet cavity oscillations (18 Hz) occurred when the Impeller 
was operating at an NPSI1 above that at which head breakdown begins, while the 
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high frequency rear chnnnol oRclllatlons (42 to 175 II«) occurred below the bond 
breakdown point. 

Badowflkl (10) discussed an Instability occurring In fou r~b laded Inducera. 

Ho found that at roduced flowratoa and a certain value of NPSH the four -’blade 
cavitation pattorn makoa a Ruddon transition to cavitation only on alternate bladea. 
Additional reduction of NPSH cauROR cnvltlos to roofltabllRh on the two noneavltntlng 
blados. An InRtablllty occurs whon tho roofltabllRhod cavltlcR roach one blado 
spacing In length. No description of a cavity Is given although roforonco to a 
aurfaco cavity model Is made. Ho states that two Instabilities have boon experi- 
enced whon the two reestablished cavities are unequal In length, ono Instability 
occurring onch tlmo a cavity grows to a length of one blado spacing. By reducing 
NPSH so that the cavltlos becomo longor than ono blado spacing stable flow was 
generally rostored. The Instability Is described as "low frequency, high amplitude 
vibration, " with no numerical valuos quotod. No mention Is mado of a propagating 
pattern. The Instability Is reported as occurring In a band of NPSH Immediately 
prior to head breakdown, with the flow stable In the head breakdown region. 

Badowskl attributes the Instability to the existence of "backflow-induced pre- 
rotation" In the Inlet llnfc. 

Etter (11) reported cavitation Instabilities on two- and three-bladed In- 
ducers designed to operate with blade cavities longer than chord length (super- 
cavltattng). Only the supercavltatlng region was Investigated. Instabilities were 
found when the two-bladed model was operated at flow coefficients below 0.090 
and the cavity lengths became less than about 240% of chord. Instabilities were 
found when the three-bladed model was operated at flow coefficients below 0. 075 
and the cavity lengths became less than about 180% of chord. In both cases the 
cavities extended back to or silently past the leading edge of adjacent blades 
when the Instabilities occurred. Flow was stable for shorter cavities. Each 
blade cavity consisted of a suction surface cavity springing from the leading 
edge and a tip vortex cavity said to be of equal length, which sometimes merged 
with the leading edge cavity. Observed Instabilities were described as "large 
magnitude, low frequency" oscillations In cavity length, with no numerical values 
given. It was noted that cavities on the two-bladed model were of equal length, 
while cavities on the three-bladed model were of unequal lengths, by as much 
as 20% but usually less than 10%, over a large operating region. In most instances 
two cavities were of one length and the third shorter or longer. In a few Instances 
all three were different. The pattern could not be predicted to occur on any par- 
ticular blades, and It was concluded that the phenomenon was not caused by geometry 
deviations among the blades. The region of observed Instabilities seemed to have 
no particular relationship to the region of unequal cavity lengths. It Is not stated 
whether tho cavities oscillated In unison. The Instabilities occurred when the 
Inducers were operating In the head breakdown region. 

Hartmann and Soltis (12) made observations of cavitation In an axial flow 
pump with 19 blades. They found that In addition to blade surface cavitation a 
largo part of the total vapor cavity resulted from tip vortex cavitation. At low 
Inlet pressure and low flow coefficients both vapor formations were observed to 
be unstable, with the tip vortex cavity fluctuating between a position close to the 
blade and a considerable distance away. At still lower flowrates the cavitation 
moved "In and out of the blade passage" such that at one Instance the passago was 
almost free of cavitation and at another almost filled with vapor. No numerical 
values are given and no mention Is made whether all blade passages were acting 
in unison. 
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Hildebrand (Id) reported cnvitntlon-inducod oscillation of the .12 engine 
oxidizer pump during engine firing testa. Frequency , aa determined by static 
proBBure probes, varied from 14 to 24 Hz depending upon onglno operating point 
nnd NPSI1. Pump Inlet and dlflehnrge proBBuroa roachod poak-to-peak ampUtudoa 
of approximately 14 N/cm 2 (20 pal) while engine chnmbor proBBure oBclllated 
with a ponk-to-ponk amplltudo of as much as 7 N/cm 2 (10 pal). It la atatod that 
ns pump Inlot proaauro was lowered n ' •* oaclllatlona, which began before tho 
head breakdown point waa reached, not. eoaao but tho frequency docroaaed. 
Addition of n hollum-flllod accumulator tionr pump Inlet roducod both the frequency 
nnd amplitude of oBclllntlon. No observations of cavity formations woro mado. 
Previous teats of tho same pump, roportod In (14), had shown that tho oscillations 
would coaso If Inlot prossuro were lowered sufficiently. It Is not stated In (Id) 
whether Inlet pressure was loworod to the samo lcvols as In the earlier tests. 

Soltis (15) Investigated cavitation formations on two different threo-blndod 
helical Inducers and an axial flow pump with 19 blades. Both a surface cavity 
with a rather distinct closure point nnd a tip vortex cavity formed. It Is stated 
that at low flowrates nnd low Inlet pressures cavitation of a highly unsteady 
nature Is observed. Significant pressure oscillations occur wltlilh the system 
and are observed as cavitation pulsations within the blado passages. In general, 
this unsteady cavitation occurs when the cavity closure point reaches Inside the 
passage formed with an adjacent blade. " The tip vortex cavitation and the surface 
cavitation both followed the same pattern of fluctuation, which " — Is evidenced 
as a chugging movement of the cavity both in and out of the blade passage. " This 
seems to Imply that the cavities on all blades oscillated in unison. However, In 
references to a film strip of cavitation on an (unidentified) helical Inducer it Is 
stated, "Here, large amounts of vapor fill some of the blade passages, while the 
other passages are almost entirely vapor free. Closer analysis of the individual 
frames Indicates that the full and empty blade passages occur in a regular pattern 
such that the cavltatlng zones rotate around the rotor at an angular speed lower 
than that of the rotor. " And in reference to the axial flow pump (19 blades), 

" — the vaporous regions have moved out ahead of tho rotor and, again, flow 
conditions are such that these regions rotate around the annulus at a speed lower 
than the blade speed of the pump. A variation of flow from this operating point 
affects the speed of rotation of the cavltatlng zone about the annulus. " No informa- 
tion on oscillation frequency or amplitude is presented. The Instabilities are 
depicted as beginning before head breakdown begins. It is not stated whether 
stable operation Is restored after head breakdown begins. 

Miller and Gross (16 ) noted cavitatlon-lnduced oscillations during testing 
of a shrouded, hubless, helical Inducer close-coupled to a centrifugal Impeller 
In water and li quid nitrogen. "High-amplitude, low-frequency Inlet and discharge 
pressure oscillations characteristic of cavltatlng pump Inducer systems were 
observed at all lower than design flowrates. At design flowrates and higher, the 
low-frequency oscillations wero observable; however, the amplitude was severely 
diminished. " Oscillation frequency and amplitude were dependent upon inlet 
pressure and flowrate, with frequency varying from 4 to 12 Hz. Inducer Inlet 
pressure reached a peak-to-peak amplitude of about 24 N/cm 2 (35 psl). Both 
frequency and amplitude decreased as inlet pressure was lowered. No 
observations of cavity motion were made. Tho data presented Indicate that the 
Instabilities occurred prior to head breakdown. 
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2.2 MODFUNG ATTEMPTS 


The reported models, employed to explain onvltntlon-lnduord oscillation In 
Inducers, fall Into the categories of rotating stall and system modeling. One 
attempt at employing n rotating atoll analyala and throe nttempln at employing 
system modollng were made. None were satla factor lly correlated with test re- 
sults. One attempt to modol the Instability In a cascade was made but agree- 
ment with test data was not achlovod, 

In a rotating stall analysis, originally developed for axial flow compressors, 
the only geometry considered Is that of the Inducer blading, with the Inlet and 
discharge piping systems assumed to have no Influence. Presence of a vapor 
cavity Is not directly considered. The major feature of an Instability Is assumed 
to be a circumferential blade-to-bladc variation or distortion In pressure and 
flowrate. '.Stter (11) used an analysis by Yeh (17) In an attempt to explain the In- 
stabilities encountered on the two-bladed and threc-bladcd superca vita ting In- 
ducers. It was found that the predicted occurrence of an Instability did not corre- 
late, even trendwlse, with the observed unstablo operating regions. Etter con- 
cluded that It was unlikely that rotating stall was responsible for the observed 
"large magnitude low frequency Instabilities. " 

In a system modeling approach all geometry Is considered; the Inducer 
blading as well as that of the Inlet and discharge piping systems. Each element 
In the system Is modeled Individually and the equations are combined Into a set 
of equations representing overall system motion. Stability Is determined either 
from linear stability analysis techniques or by simulation on an analog computer, 
depending upon the form of the equations. The presence of blade cavities Is 
directly considered through application of the law of continuity of mass to the 
fluid entering and discharging from the Inducer. It Is assumed that discharge 
flowrate can be different from Inlet flowrate because total vapor cavity size can 
change. The major feature of an Instability Is assumed to be oscillation In Inlet 
and discharge flowrates, with no variation In pressure and flowrate between any 
two blade passages. This Implies that the cavities are equal In size and move In 
unison. The three reported modeling attempts have all neglected tip vortex cavi- 
tation, treating only the blade suction surface cavity. 

Gross (18) modeled a shrouded, hubless, helical Inducer close-coupled to a 
centrifugal stage operating In a water test facility. The facility exhibited low fre- 
quency (3 to 14 Hz) high amplitude (up to ±70%) Inlet pressure oscillations over a 
wide range of Inlet pressures and flowrates. Only the region above the head 
breakdown point was Investigated. Both frequency and amplitude decreased as 
Inlet pressure was lowered. A system model of the facility was programmed on 
an analog computer. Calculation of cavity volume was based on an assumed 
parabolic profile with the height computed from the Strlpllne and Acosta (19) two- 
dimensional cavity model and the length-to-helght ratio baaed on an experimental 
point. The resulting model would not exhibit an Instability. 

Sack and Nottage (20) modeled a four-bladed helical Inducer operating In a 
liquid oxygen test loop. The loop Is stated to have been unstable, without speci- 
fying the exact operating region of Instabilities. A typical frequency mentioned 
Is 5 Hz, and discharge pressure oscillation Is shown to have reached an ampli- 
tude of ±8 N/cm^ (±12 psl), A system model of the loop was programmed on an 
analog computer. Calculation of cavity volume was based on an assumed circular 


arc profllo tangent to tho fluid Incldonoo anglo, with tho height computed from 
the Stripling and Acosta cavity modol (19), Operating regions both abovo and 
bolow the head breakdown point were Investigated with the model, The model 
oxhlbited Instabilities in both regions for certain l)andH of Inlet pressure, Tho 
only correlation with test results presented are the statements that "frequencies 
and wave shapes were representative" and, "the real system ns tested was ob- 
served to limit cycle when unstable. The simulation was noted to have a limit— 
cycling amplitude which also followed observed trends with cavitation number. 

The simulation had a limit-cycle amplitude higher than the lest inducer and was 
unstable for a wider range of cavitation number." 

A three-bladed helical Inducer operating In a water loop was modeled In 
(21) on an analog computer. Calculation of cavity volume was based on the exact 
profllo computed by tho Stripling and Acosta cavity modol (19) , It Is stated that 
the rosultlng system modol was unstable, but work was terminated and no correla- 
tion with test results was presented. 

Acosta (0) modeled a cascade in which unisonous cavity oscillations 
occurred. The details of the model were not presented, but it apparently was a 
system model. It is stated, "It is easy to show that the dynamic system consist- 
ing of the nozzle flow, the cavilating flow in the working section, plus some 
allowance for the inertia of the diffuser flow leads to a third order system that is 
inherently stable provided only that the total pressure loss due to cavitation In the 
working section increases as the ambient pressure decreases. Total pressure 
loss across the cascade was computed as a function of time from a knowledge of 
the instantaneous velocities and pressure levels upstream and downstream of 
the cascade. In this calculation the effect of the accelerating main stream from 
the point of measurement of the pressures was taken into account. This is a diffi- 
cult calculation to make because of the spikiness of the velocity traces; these 
in turn caused large changes in pressure between the cascade and the reference 
point due to the large acceleration of the main stream. Nevertheless, it does 
appear that the total pressure loss Is least when static pressure upstream of the 
cascade is greatest, . , " 

2. 3 DISCUSSION OF PREVIOUS RESULTS 

Experimental observations of cavitation-induced oscillation reported in the 
literature do not present a completely clear picture of the phenomenon. The ob- 
servations on individual hydrofoils (2, 3, 4) and cascades (5, 6) strongly Indicate 
that the blade surface cavity is a majorfactor. However, Wood (9) and Hartmann 
and Soltis (12) indicate that additional cavity formations in the tip vortex or insti- 
gated by the tip vortex may also be involved. The hydrofoil and cascade tests 
indicate that unisonous cavity oscillations occur. However, some of the inducer 
tests show circumferentially propagating cavitation (8, 15), or alternate blade 
cavity oscillations (8, 10) . 

There could easily be more than one vapor cavity involved in the observed 
instabilities. In general, in both an inducer and a cascade, there will exist both 
a blade surface cavity springing from the leading edge and a tip vortex cavity. 

The tip vortex cavity generally forms at higher NPSH than the surface cavity. It 
is reasonable to assume that as NPSH is lowered the tip vortex cavity could form 
and begin to fluctuate in position (or form and detach) while the blade surface 
cavity remained steady or had not yet formed. This could occur because of the 
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the mechanism of the phenomenon It was found not to favor any certain bladea, 

If cavitation-induced oscillation oceura with the cavltlea moving In unlnon during 
alternate blade cavitation, Ihen aa the cavltlea collopae during one cycle of oaollln- 
tlon the tendency to re-form on the other two bladea may be greater than the tend- 
ency to re-form on the an me two bladea. The combination of cavitation-induced 
oaclllatlon and alternate blade oavllatlon could then appear aa n circumferentially 
propagating cavitation pattern. 

From consideration of the obaorved cavitation phenomena, It in concluded 
that modeling of cavitation-induced oscillation should consider as major features 
(1) a suction surface cavity starting at the leading edge of each blade, and (2) 
oscillation of these cavities in unison. Tip vortex cavitation caused by leakage 
flow through the blade tip clonranco space should bo considered of secondary 
Importance in Initial modeling. Propagating cavitation should be considered a 
different phenomenon from cavltatlon-lnduced oscillation. 

Experimental evidence to support the above conclusions is scant, primarily 
because of the lack of a systematic investigation of cavitation formations and their 
relationship to Instabilities. Oscillation of the blade surface cavity is the major 
feature mentioned in the individual hydrofoil tests (2, 3, 4), and in the cascade 
tests (5). In the inducer tests a surface cavity is either stated to be present 
(9, 11, 12, 15), implied to be present (8, 9, 10) . or no visual observation was 
made (13, 16). There is no reported case of an instability in which it is known 
that there were no blade surface cavities. 

The experimental evidence indicates that although there may be more than 
one mode of cavity motion - propagation as well as oscillation - the most likely 
mode occurring during most Instances of low frequency, high amplitude cavitation- 
induced oscillation is that of nonpropagating oscillation with all cavities moving 
in unison. This is inferred from the individual hydrofoils tests in which no propa- 
gation is possible and from the absence of a propagation pattern in the cascade 
tests. The strongest evidence that propagation is not the major feature of the 
instability in question is the overall flowrate oscillation indicated by chamber 
pressure oscillations in (13). A purely propagating cavitation pattern would not 
be expected to cause an oscillation in overall inducer flowrate. 

It is likely that tip vortex cavitation can influence stability and under certain 
conditions become an important factor. The presence of tip vortex cavitation 
during an instability was noted in the cascade tests (5, 11) and most of the inducer 
tests (13, 9, 11, 12). However, the tip vortex cavity "was treated as a secondary 
consideration in thiB program because there exists no verified model for computing 
vapor volume of the tip vortex cavity. The only known published model is by 
Ghahrer.ianl (22). The influence of the tip vortex cavity con be estimated by the 
success In predicting instabilities when It is Ignored, allowing an estimate of the 
level of effort which should be dlrocted toward a tip vortex model. 


SECTION a 

APPROACH TO INSTAniUTY ANAT.YflIH 


After consideration of the reported descriptions of cavity motion during 
cavitation- Induced osclllntlon In Inducer, ciiscndo, mid hydrofoil tests, and of 
previous modeling attempts, wo concluded that tho ay atom modeling approach 
offered the grealont potential for the prediction of Instabilities. There In no 
evidence tliaf rotating stall or hlado-to-hlade dlatortlon (which are not con- 
sidered In a system model) wore significant factors In most known Instances 
of cavitation- induced oscillation. The system modeling approach can, however, 
qualitatively < xplaln all reported Instances of Instability and appears to bo 
representative of the actual flow process. 

:i. 1 MECHANISMS THAT CAUSE IN'STAIUUTIES 

Two basic mechanisms can create on Instability under cavltatlng condi- 
tions: 

1. Tho continuity mechanism, which results from tho relation- 
ship of cavity volume to Inducer Inlet flowrate, and requires 
only tho presence of cavitation cavities on tho Inducer blades 
for an Instability to be possible. 

2. The performance mechanism, which results from tho relation- 
ship of Inducer headrlse to Inlet flowrate, and requires that 
cavitation affect the inducer performance for an Instability 

to be possible. 

The first mechanism begins at a cavitation number well above head break- 
down, and can be physically understood by considering that the blade suction 
surface cavity streamline constitutes a flow boundary. Since the streamline 
can change position In response to flowrate changes, It can be considered a 
movable wall. A momentary Increase In flowrate reduces the fluid Incidence 
angle, causing a reduction In cavity volume and a corresponding Increase In 
the local space available for liquid. This Is Illustrated In figure 1. Continuity 
requires that Inlet flowrate nccelerate to fill the newly created ttquld space, 
but this reduces Incidence angle and causes a further collapse of the cavity. 

Inertia of the fluid mass In the Inlet line causes the accelerating flowrate to be 
accompanied by a reduction In static pressuro at the Inducer Inlet. The reduced 
Inlet pressure acts to counter cavity collapse, since lower static pressure causes 
a cavity volume to Increase, and at some point the static pressure reduction Is 
sufficient to prevent further cavity collapse. Fluid Inertia then causes an over- 
shoot to a reverse trend, during which tho Increase In flowrate (with accompanying 
decrease In pressure) results In an Increase In cavity volume. The flowrate 
Increase Is eventually halted an i a deceleration begins with an accompanying 
Increase In pressure. The deceleration continues back through the original 
operating point, where a decrease In flowrate causes an Incroaso In cavity 
volume, and at some point the accompanying pressure Increase causes tho 
cavity to cease growing and begin to collapse again. Inertia causes an under- 
shoot, the deceleration Is halted, and the flowrate acceleration boglns again. 

The path of oscillation Is depleted In figure 2. 
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CAVITY VOLUME 




Thlfl In descriptive of tho moohnnlflm that cauHon Inducor Inniribl litlon 
at cnvltntlon numbers above tho bond breakdown point* Tho mechanism oIro 
exists nt on vl lotion numborn bolov' tho point nt which howl breakdown begins, 
but howl breakdown Introduces n stabilizing offoot. In <ho howl bronkdown 
vpglnn tho lowering Inlot proRnuro neoompnnylng u flownito acceleration not 
only note to counter onvlty oollnpno, but iilno lo lower Inducor dlsehnrge I'l’en- 
ru re (through pump gain), which tondn to roduoo tho flownito thiil oim bo forced 
through the discharge system. Tho deceleration of dlsehnrge llowrnlo helps to 
fill tho nowly crontod fluid wpiico and reduces tho accolonillon rouulrodof Inlot 
flowrate, 'I'll In explains why tho osellliitloim have gonorally boon noted to 
dlmlnlRh once howl breakdown begins. 

Tho second mechanism begins In tho bond breakdown region and Involves 
tho way cavitation can affect Inducor headrlne* If an Inducer In operating nt a 
point at which a momentary IncronHo In flowrate onuses nn Increaso In hoadrlac, 
as shown In figure 3, tho Instantaneously higher head nvallablo becomes greater 
than the head required to maintain flow and tho flowrate Increases further. 
Accompanying tho flowrate Increase, Increasing frictional resistance In tho 
Inlet line lowers Inducor Inlet static pressure, which tends to Increase tho 
degree of head breakdown, and the Inducor operating point shifts to lower values 
of headrlsc on the head vs flow map. Doth offocts tend to reduce the hoad 
available. The Increasing flowrate also causes greater pressure drops across 
discharge restrictions, Increasing tho head required to maintain flow. Those 
effects gradually stop the flowrate Increase. An inertial overshoot causes 
head available to drop below hoad required to maintain flow, and flowrate begins 
to decrease. At this point It Is speculated that hysteresis In tho Inducer system 
keeps head available below head required until ilowrate drops below the original 
operating point. An Inertial undershoot then causes head available to Jump 
above head required, and flowrate begins to Increase again. The path of 
oscillation Is depicted In figure 4. 

Both destabilizing mechanisms are influenced by the hydraui.. * system In 
which the Inducer Is operating, because tho system haB an effect on the exact 
relationship between pressure and flowrate. The system can therefore Influence 
the point at which oscillations begin and the frequency and amplitude of oscilla- 
tion. For this reason, the entire hydraulic system must be considered in an 
analysis of cavltatlon-lnduced Instabilities, 

Both mechanisms can occur In Inducers, while only the latter Is believed 
to occur in a cascade. In an Inducer, an Increase In flow decreases Incldenco 
and cavitation No. Tho Incidence decrease generally Is a stronger offcct on 
cavity size than the cavitation No. decrease, so that cavity size decreases. In 
a cascade, Incidence Is fixed so that on Increase In flow causes an Increase In 
cavity size. Thus the continuity mechanism would not be expected to cause an In- 
stability In a cascade. The performance mechanism can causo u cascade Instability, 
As cavitation number Is lowered, the experimental observations show that steady- 
state lift and drag forces on the cascade remain constant until a point Is reached at 
which both forces begin to change. Further lowering of cavitation number causes 
first one force and then the other to usually Increaso, roach a peak, and then 
rapidly decrease (2, 6, 6). The reported cascade (and hydrofoil) Instabilities 
occurred In this lift and drag breakdown region. Since the lift and drag forces 
are forces exerted by the cascade on tho fluid, they create pressure-area forces 


fit tho nan carlo Inlet anil exit planon, The prosnuro-aroa forces In turn not on 
tlio flow n pp roa olilng' nnd leaving tho on sonde, If an Inoronnn In flow oauaoa tho 
lift, drug, nnd pronfiuro-iiron forooa to ohango no an to nooolorato tho flow, 
tho moulting flow Increase will further ohnngn tho forces nnd nn unntnblo trend 
dovolopn, Thin trend will continue until n change In flow or pronauro oaunen the 
force ohnngen to rovorno dlrootlon. 
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Figure 3. Hoadrlse Response to Flowrate; Perform- 
ance Destabilizing Mechanism Below Head 
Breakdown 


FD 82373 




Hcudme Required 
Headrise Available 


FLOWRATE 


Figure 4. Path of "Performance" Mechanism Instability FD 62374 


It should bo pointed out that those rooohnnlsmn are not modeled Individually. 
Only n single model I a developed, end thnt model contains both mnehrmtamr of mi 
Instability. Identifying the sopnrfttn mechanisms wnn on nlo In determining the 
offsets required In n model. 


a, 8 DKFINITinN OF llWjUIHKn KXl'KIUMMNTAI. DATA 


There' nre two functions (iipnrt from the readily pre Itelahle Inlei and 
discharge line ehn meter! iitlon) that require aeeurate quanlltatlve definition 
for Hie analytle prediction of the oeeurrenoe and eharaeterlnllen of aelf-lnduneil 
Instabilities In eav Hating Induoora, f rhe flrat of theme In eavltatlon eavlty 
volume no a function of Inductor Inlet pressure nnd flowrate, and the nooonil In 
Indueer head rlno an a function of Inlctt pressure and flowrate. Thene functions 
are associated with the "continuity” and "performance" Instability mechanisms, 
respectively. 

At present, the cavity volume function la considered the most Important 
beenuBO It affectB tlu* occurrence of InaUibllltloo In the usual Inducer operating 
regime - reasonably above the head fa lloff pressure. It should bo noted that 
total cavity volume (blade suction surface and tip vortex cavitation) must bo de- 
fined In a rigorous analysis. Tlio fact that the continuity Instability mechanism 
requires that cavity volume bo sensitive to flow race, and that a qualitative con- 
sideration of tip clearance cavitation Indicates that Its volume may bo much less 
sensitive to flowrate than Is suction surface cavitation volume, led to the con- 
clusion that an Initial analytic modeling attempt would provide a reasonable degree 
of accuracy if only the suction surface envlty wore considered. 

Suction surface cavity volume for arbitrary inducer goometrles and flow 
conditions Is predictable with the computer program of (1), but the predictions 
did not provide the accuracy roqulred because of the poor reported correlation 
of predicted with measured cavity length data. Additional experimental data 
wore therefore required on the geometry of blade suction surface cavitation as 
a function of Inlet pressure and incidence angle (flow coefficient) for refinement 
of the cavity model from (1). 

A second set of experimental data were required to demonstrate the adequacy 
of the analytic instability modeling technique, to provide a qualitative indication 
of the accuracy of the suction surface cavity model, and to provide a r« allta- 
ttve Indication of the relative importance of Up vortex cavitation. 


The experiments that were planned and conducted to satisfy these require- 
ments are reported in the next section of this report. 
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Two experiments woro conducted: n cascade experiment and an Inducer 
experiment. The primary objective of the cascade experiment wns to provldo 
cavity goomotry data, and a secondary objective was to provldo Instability data 
Cur correlation with analytic predictions. The objective of the Inducer oxporl- 
mont was to provide Instability data for correlation with predictions. The experi- 
ments and their results aro discussed In this Section, 

4.1 CASCADN RXPKIUMRNT 

4.1.1 Facility 

The cascade experiment was conducted In the water cascade tunnel at The 
United Aircraft Rosoarch Laboratories. A drawing of the tunnel Is shown In 
figure 8. Tho tunnel was originally designed and constructed to permit the 
measurement of two dimensional performance and cavitation Inception points of 
various cascade configurations. These previous programs are reported In (7). 

The tunnel Is vertical and Is oriented so that the Inlet plane of the cascade Is 
" horizontal, which eliminates gravitational hydrostatic pressure gradients along 

the length of the cascade. 

Water Is circulated In the closed-loop tunnel by four pumps, which are 

V located below floor level to Increase available pump Inlet head. The flow Is dls- 

V charged from the pumps, diffused and turned In a system of ducts, and settled 

In a rectangular chamber containing both a honeycomb flow stralghtener and 
graded-poros Ity screens for reducing large-scale turbulence, A subsequent 

n guide vane section aligns the flow and delivers It to the cascade Inlet nozzle. 

The flow Is accelerated by the Inlet nozzle, then passed through the cascade 
test section, and finally discharged Into the plenum tank from which the flow Is 
drawn Into the circulating pumps. 

Four Interchangeable cascade Inlet nozzles are available for Inlet flow 
angles of 0,87, 1,06, 1,22, and 1,31 rad (60, 60, 70, and 76 deg', as measured 
rS from the nozzle axis to a line that Is normal to the cascade Inlet plane. The 

1.31-rad (75-deg) nozzle was used In this program because It most nearly approxi- 
mates the blade angles typical of Inducers. Installation of the 1.31-rad (76-deg) 
Inlet nozzle with the test section assembly Is shown In figure 6. A section In the 
top wall of the nozzle Is transparent to permit observation of cascade Inlet flow 
conditions from a position above the test section. The cascade section Is normally 
assembled with porous side walls In place of the Illustrated transparent acrylic 
walls to provide for removal of the cascade boundary layer. The transparent 
— walls wore used In this program to permit visualization of cavitation profiles. 

The resulting partial loss of cascade two dimensionality was not considered signifi- 
cant to this program In view of other experience with oavltatlng cascades (2, 4, 

6, 6). 


c 


1_L 


— 


1? 


i 





1 


J, 


o 



Figure 5. Water Cascade Tunnel FD S3038 



Figure 6. Hydrofoil Cascade Installed In Water FE 113784 

Tunnel 

A drawing of the cascade test section area Is shown In figure 7. The test 
section Inlet Is rectangular, having a 15.24-cm (6-ln.) span and a 91.44-cm 
(36-ln.) length. The test section was designed to establish both a uniform flow 
along the cascade Inlet plane and a periodic (blade-to-blade) flow downstream of 
the cascade. These are necessary conditions to ensure that the flow about each 
blade Is Identical; accomplishment of these goals provides a flow that Is representa- 
tive of the flow through an Infinite cascade. Approximations of flow uniformity and 
periodicity are achieved by removal of the wall boundary layers and by contouring 
the ca^ade end walls to minimize their effect on the mainstream flow. The 
boundary layers which develop along the walls of the Inlet nozzle are removed 
upstream of the test section by means of step-type slots located on all four walls 
of the nozzle. Flow Intercepted by each of the four step slots Is ducted through 
Individual throttling valves to one of the four pumps. 

Adjustments of flow streamlines (to obtain uniformity) are accomplished 
by means of variable geometry end walls located at each end of the cascade. 

These end walls are comprised of throe sections; a flexible wall that connects 
the rigid Inlet nozzle to an adjustable end wall, an adjustable end wall that Is 
analogous to one surface of a cascade blade, and a tailboard that extends down- 
stream from the adjustable end wall. Actuation of these end-wall sections enables 
Independent adjustment of (1) the gap between the end blades of the cascade and the 
adjustable end wall, (2) the angle of the adjustable end wall, and (3) the angle of 
the tailboard. The convex, flexible end wall Is porous to permit removal of the end 
wall boundary layer, thereby decreasing any tendoncy toward separation from this 
surface. 









The cascado tost section naaombly Ifl flupportod from the fixed Aide of the 
3, 06-meter (10— ft) diameter plonum tank and completely encloflod within tho tank. 
Tank wator level during testa Is Just above tho teat flection. Tho tnnk hna an 
operating pressure capability of 0.7 to 09 N/cm 2 (1 to 100 pain), Tho throo main 
flow pumps and ono boundary layor pump arc fabricated from zlnc-froo bronze 
and stainless steel. Each pump Is driven by a 7457-watt (10-hp) motor and Is 
capablo of delivering water flows at a rate of 6.44 m 3 /mln (1700 gpm) with a 
head rise of 4.88 m (16 ft). 

Contamination of the water used In tho facility Is minimized by special 
water processing equipment. In addition, tho facility is constructed primarily 
of stainless steel to avoid solid particle shedding. Solid particles contained In 
the water aro romoved by various filters, one of which provides continuous 
3 x 10" 6 m (Sfx) filtration at a flowrate of 0.38 m 3 /mln (100 gpm). Dissolved 
minerals are removed from the test water by a commercial lon-exchange-type 
demtnerallzer that provides water comparable to distilled wator In conductivity, 
a measure of dissolved mineral content. The conductivity of the water from this 
unit was approximately 1 micro ohm/cm. Dissolved gases may be removed by 
a cold-water deaerator that can reduce the gas content to 3 ppm. The deaerator 
was not used In this program, however, as the natural deaeration resulting from 
the low test pressure was found to be more efficient. 

4.1.2 Cascade Configuration 

The test cascade was designed to be as nearly representative of Inducer 
blading as practical within the limits of existing tunnel hardware. The selected 
cascade features are listed In table 1 and a drawing of the cascade Is shown in 
figure 8. The flat plate profile Is the cascade equivalent of helical Inducer blades. 

A cascade solidity of 2.0 was selected as being representative of Inducers (but on 
the low end of the usual range) to permit large, easily measurable cavities to 
form. Hydrofoil chord was established near the upper limit of the test section 
capability so that the selected solidity would result In as large a blade spactng, 
and as high a cavitation cavity (for measurement), as possible. Leading edge 
wedge angle was determined by the requirement that the blade suction surface 
clear the free uow streamline at minimum Incidence, and the leading edge was 
made relatively sharp (Inducer design guide lines for best suction performance). 

The trailing edge configuration was unimportant In this program, where cavitation 
and not performance was to be studted, so the trailing edge was merely rounded 
to simplify manufacturing. Hydrofoil thickness was selected to make blade bowing 
negligible, yet not so thick as to result In unrepresentative blade blockage effects. 
Hydrofoil span was set solely by the test section capability and the number of hydro- 
foils resulted from the selected soltdlty, chord, and the length of the cascade test 
section (91.4 cm, 36.0 In.). Inlet flow angle was set by the existing tunnel 
facility nozzle. The nozzle with the highest angle (1.31 rad, 75 deg) was used 
to approximate, as nearly as possible, the blade angles of a typical Inducer 
(1.36 to 1.48 rad, 78 to 85 deg from axial). 


Table 1, Cascade Foaturoa 


Hydrofoil Profile 
Solidity 

Hydrofoil Chord 
Hydrofoil Spacing 
Loading Kdge Shape 

Loading Kdge Radius 
Trailing Kdge Radius 
Hydrofoil Thickness 
Hydrofoil Span 
Number of Hydrofoils 
Inlet Flow Angle 


Flat Plato 

2.0 

IB. 2 om (6.0 In.) 

7.02 cm (3.00 In.) 

0.09-rad (5-dog) wodgo angle cut 
on suction surface 
0.025 cm (0.010 In.) 

0,25 cm (0.10 In.) 

0.61 cm (0.20 In.) 

15.2 cm (6.0 In.) 

11 

1,31 rad (76 deg) 



Figure 8. Cascade Configuration 


FD 62668 


The hydrofoils were machined from AlSI-type 304 stainless steel and were 
welded to studs, one at either end of the hydrofoil. A photograph of a com^eted 
hvdrofoll Is shown In figure 9. The hydrofoil surfaces were not polished after 
machining because the Tact cavitation 

of Interest. Surface roughness measurements of a typical hydrofoil wer . 

Rouehness (RMS) 


Suction Surface 
Wedge 

Flat 


3.8 gm (150 /min.) chordwlse 
0.6 /m ( 26 Mta.) spanwlse 
1.3 (60 /Uln.) 


Pressure Surface 


6.4 jum (260 /uln.) 


The cascade hydrofoils were assembled In the tunnel test section between 
acrylic plastic sidewalls, as shown In figures 10 and 11. The hydrofoil studs 
extended through holes In the sidewalls so that the foils pivoted about tae stud 
centerline when Incidence angle was changed. Foil end clearance was Just 
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sufficient to permit smooth pivot operation and was approximately 0,013 cm 
(0,000 In.) for all foils (total for both ends), A series of rbferonec linos, one 
at oach hydrofoil location, wore scribed on tho far sidewall prior to assembly 
to facilitate sotting tho foils at a uniform Incldonco angle, Tho studr were loolcod 
to an actuation llnkago with setscrews along tho scribed lines so that subsequent 
actuation of tho llnkago, which was accomplished romotoly, movod nil the foils 
simultaneously through the snmo nnglo. 

4.1,3 Instrumentation 

Cascade tunnel Instrumentation was Installed to provide data on (1) cavity 
geometry ns a function of operating condition, (2) cascade pressure loss, and 
(3) tunnel dynamic pressures. 

Cavity geometry measurements wore obtained from photographs taken 
looking through the acrylic cascade sidewall at the hydrofoil profile and through 
the window In the top of the nozzle at the hydrofoil suction surfaces. The photo- 
graphs were taken from outside tho tunnel through portholes In the plenum tank, 

A spccdgraphlc 4x5 camera using Kodak Trl X "llm and a Calumet 4x5 long 
focus view camera using Polaroid type 47 film were used with Illumination pro- 
vided by a Graphic Stroboflash - 1 strobe light. 
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Figure 9, Cascade Hydrofoil 


FE 119531 









Cascade Inlot total pressuro, static prossuro, and velocity prossuro wore 
obtained from a Kiel probo at the entrance to the Inlet nozzle and a wall atatlc 
tap approximately 4, 0 cm (1. 8 in. 1 upatroam of tho caaonde loading edge plane, 
Each parameter waa road from Ita own mercury manometer, A aerloa of wall 
atatlo tapa apaeod at 2 , B-em (1-ln, ) Intorvala along tho length of tho eaaeade 
at the eaaeade Inlot and dlachargo and a Kiel probo approximately 25 cm (10 In. ) 
downstream of tho oaacado wore oonnootod to a 204-cm (100-ln, > tall, BO-tubo 
mnnomotor to Indicate onacado flow uniformity and dlachargo totnl pressure, 
rospoetlvoly. Thoao meaauromonts wore later found to bo Inaccurate undor 
eavltntlng conditions boenuao of air bubbles that dovolopud In tho long mnnomotor 
lines. Wator tomporaturo was moasured with an Immersed morcury thormomotor. 
Plonum tank ambient pressuro (abovo tho wator froo surfaco) was moasurod with 
a pressure gage. Maximum Instrumentation error Is estimated nt +0. 085 cm ilg 
(*0. 025 In. Ilg) for the manometer measurements, il. 25 cm Ilg (i0. 5 In. Hg) for 
the pressuro gage, and h0.3°K (i 0.5°F) for the thermometer. 

Dynamic prossure data were measured by means of four Klstler model 601 H 
piezoelectric transducers with Klstler model 668 charge amplifiers. The trans- 
ducers were located In the tunnel as follows: 

1. In the large section at the nozzle entrance, 238 cm (94 In. ) 
upstream of the cascade leading edge plane 

2. In the convergent portion of the nozzle, 165 cm (65 In. ) upstream 
of the cascade leading edge plane 

3. At the cascade Inlet, 1 cm (0. 4 In. ) upstream of the leading 
edge plane 

4. At the cascade discharge, 46 cm (18 In. ) downstream of the 
leading edge. 

These locations are Indicated on figure 5. Data were recorded using a Honeywell 
Model 6600 magnetic tape recorder for the first 18 test points and a vlslcorder 
Model 1108 oscillograph for test points 19 through 41. 

4. 1. 4 Test Procedure 

The cascade was Initially tested at high (noncavltatlng) pressures and 
attempts were made to obtain uniform Inlet and exit flow distributions through con- 
trol of boundary layer flow and positioning of endwalls and tailboards. The degree 
of flow uniformity was determined by inspection of the manometer pressure pro- 
files from the Inlet and discharge static wall taps. It was necessary to adjust 
boundary removal flow, endwall setting, and tallborad setting to maintain an 
acceptable uniformity over a range of Incidence angles from 0. 09 to 0. 26 rad 
(5 to 15 deg). Under cavltatlng conditions, manometer pressure profiles could 
not be used to determine the uniformity of flow because of uncertainties In the 
Individual pressure values caused by air bubbles trapped In the long lengths of 
small diameter tubing +hat ran to the manometer board. Under these conditions, 
flow uniformity waB assumed when cavity lengths on the various foils were 
equal. 
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Endwall nnd tailboard Rottlngn that ontnbllflhod tho moat uniform flow during 
nouenvltatlng oporntlon did not noooRnnrlly maintain uniformity aftor envttntlon 
Inception. Adjuatmontfl woro ununlly required after the Rotting of envltatlng teat 
points, and the attainable degree of uniformity became progrenalvely worae aa 
Inlet pronnuro wan lowered nnd the cavities boenmo larger. 

During noneavltntlng tunnel oporntlon, miscode Inlot and dlachargo spun- 
wtfie velocity profiles wore monHurod at mid length along tho onsendo to obtain 
an Indication of cascade two-dimensional lly. The measurements were taken with 
tho hy drofoils set at D, 17-rad (10-dog) Incidence angle. ) An upstream boundary 
layer thickness of approximately 1 cm ( 0. 4 In. , 0. (17% span) nnd a downstream 
thickness of 4. 5 cm (1.75 In. , 29% span! woro measured. Thleknoss Is defined 
ns the distance to 90% of froestream voloelty. These valuos wore considered 
acceptable for the purposes of this Investigation. Tho velocity probes could 
only bo used during noneavltntlng operation (cavitation developed on the probes), 
and they were removed boforo the cavity geometry test points were run. 

Hydrofoil pressure surfaco cavitation was noted at Incidence angles from 
0. 09 to 0. 23 rad (5 to 13 deg) and to a small extent at Incidence angles of 0. 24 
and 0. 26 rad (14 and 16 dog). The cause of the pressure surface cavitation was 
not determined conclusively but It was believed to be caused by deviations of 
the Inlet flow angle from the presumed 1. 31-rad (76-deg) nozzle angle. Such 
deviations could have been causod by cavitation choking In the stop-slot boundary 
layer removal system, with consequent loss of the Infinite cascade simulation. 

Flow angle, and Incidence, would vary from blade to blade In such a situation 
and this variation would be manifested In nonuniform cavity geometries (that 
were observed) and In pressure surface cavitation. 

A tuft was 'nserted In the cascade Inlet flow at one location approximately 
1 cm (0. 4 In. ) upstream of the leading edge plane to gain a qualitative Indication 
of the magnitude of any flow angle variation. The location Is probably within 
the region of normal cascade Inlet flow adjustment, and, as Buch, was only used 
to roughly Indicate the freestream direction. The tuft was near the top of the 
acrylic sidewall and could not be located any appreciable amount furthor upstream. 

A flow angle of approximately 1. 25 rad (66 deg) was Indicated by the tuft at a 
typical cavltatlng test point, which Is 0. 16 rad (9 deg) less than the nozzle angle. 

The apparent 0. 16 rad upstream flow adjustment was somewhat larger than would 
be expected from normal two-dimensional considerations (0,03 to 0. 10 rad In 
potential flow calculations for a similar cascade), and the tuft angle, therefore, 
was considered to Indicate that flow adjustment was taking place upstream of 
the cascade. 

Modifications to the facility, which may have alleviated the upstream flow 
adjustment problem, such as redesign or elimination of the step-slot boundary 
layer control system, were beyond the scope of this program and could not be 
accomplished, Alternatively, the cascade was operated at an Indicated Incidence 
angle that was sufficiently high to produce suction surface cavitation exclusively, 
so that the resulting test data would give a qualitative Indication of the required 
cavity geometry vs Incidence angle and Inlet pressure functions. Step-slot boundary 
layer (b/1) flow was removed at full capacity for all test points because of previous 
experience which Indicated that excess b/1 flow had a lesser effect on two- 
dlmenslonallty than insufficient flow. In addition, several points were taken with 
no b/1 removal flow to permit an evaluation of the effects of b/1 bleed. Additional 
Inducer test data wore available from (1) to permit quantitative evaluation of the 
data. 
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CriRtmtlo clfitn wore taken nt Indicated Incidence angles of (), 28 , n, ai, and 
0. 3D rad (10, IS, and 20 dog). Theao worn tho lownnt valuoa at which no pressure 
Hurfnoo cavitation occurred and which produced reasonably a l/, ad miction surface 
cavities, All throe Incidence angina were looted with an Inlet velocity of approx- 
imately 0,130 m/p (21,0 ft/nne), and the 0, 2H»rnd Incidence wan repeated with a 
velocity of approximately 13, (M m/a (18,0 fl/seo). Blade chord Hnynoldn numhera 
wore approximately 1,<l x 10 ;i and 1.2 x 10”, reapectlvely, At each Incidence 
angle, Inlet pressure woh reduced to a level whore the caaoade waa supor- 
onvltntlng (cavity length waa one chord or longer), and preaaure waa Incroaaod 
In approximately ton atepa until the cavitation dlaappeared. At each Inlet preaaure 
teat point endwall and tall aottlnga were adjuated until blado-to-hlmle cavitation 
waa na uniform nR possible, photographa of the enacado profile and auction aide 
woro taken, n tape or an oaelllograph recording of the dynamic preaaure trnnaduoor 
algnals was made, operating condition dutu woro manually recorded, and visual 
observations concerning tho cavitation were noted, 

Tho water was deionized before testing was started In an lon-oxchnnge-typo 
demlnorallzer and was continuously filtered nt a rate of 0.38 cubic meters por 
minute (100 gpm). A low all* content was achieved by maintaining tho system 
nt pressures below 5 centimeters (2 In.) of mercury absolute at all tlmos excopt 
during the actual testing periods. Air content was not measured, but this pro- 
cedure rosultod In values of loss than 3 parts por million In provlous progrnms. 

4. 1. 5 Data Reduction 

All measurements and data reduction calculations were done In the English 
system of units (pressure in Inches of mercury, lengths In Inches or feet, etc.), 
and converted to the SI for presentation in this report. The originally measured 
or calculated data are shown as the secondary units in tables and figures. Sym- 
bols are defined in the Appendix. 

Cascade inlet velocity (Vj) was calculated from the measured difference be- 
tween total pressure upstream of the inlet nozzle and static pressure at the 
cascade inlet plane (q^) 



Tho cavitation number (k) was calculated from the cascade inlet static pres- 
sure (pj), the vapor pressure of water (p v ) at the moqsured temperature, and the 
inlet djmamic pressure (q^) 



Cascade pressure Iobb (D) was defined as the difference between total pres- 
sure at the entrance to the inlot nozzle and plenum tank static pressuro and was 
calculated from the measured data as follows: 



An uncertainty estimate was performed for the cascade pressure drop and 
cavitation number calculations since inflections In their relationship were slgnlf- 
leant to tho Instability correlation, Tho estimated measurement error In each 
of the values used In the calculations wasi 


Max Error, la, H g 


Inlet Static Pressure, Pj 

i 0.03 

Fluid Vapor Pressure, p v 

• 0.02 

Inlet Velocity Pressure, Hj 

i 0,03 

Tank Static Pressure, p^ 

i 0,00 

Atmospheric Pressure, p ft 

Negligible 


Tho uncertainty analysis was performed ns detailed in (23) for nominal data values 
that '..ore representative of data points over tho tested range of cavitation number. 
The maximum uncertainty war, i 2% for cavitation number (hi and i 9% for pres- 
sure loss. 

Cavitation cavity length and height were measured from tho cascade suction 
surface and profile photographs, respectively. Scalo for the profile photographs 
was obtained from a grid on the transparent cascade sidewall and suction surface 
photograph scalo was obtained from the knowa length of tho hydrofoil leading edge 
wedge. Both sets of measurements were corrected for viewing anglo. Height 
measurements were taken from the suction surface in a direction notmal to the 
surface, and length measurement, were taken parallol to the suction surfaco 
Both measurements were divided by foil tangential spacing (V. 6 cm, 3.0 in. ) for 
correlation with predictions. 

Dynamic pressure amplitude and frequency were read from an oscillograph 
for the nozzle, cascade inlet, and cascade discharge tranducers. The reported 
amplitude is the measured peak-to-peak amplitude of the total signal, which was 
not a simple sine wave, and frequency was calculated from the period of the most 
prominent frequency. The reported frequencies are approximate because the 
signal was erratic. 

4,1.6 Test Results 

A list of test points showing the controlled parameters, fluid conditions, 
plenum tank static pressure, calculated pressure loss, vapor pressure, Inlet 
velocity, and cavitation number Is shown In table 2. Forty-oue test points were 
run at the listed Incidence angles, Inlet velocities, and Inlet static pressures. 

All points wero run with the tunnel boundary layer bleed slots In operation at 
full capacity, with tho exceptions of points 19, 31, and 41., Those were run with 
no bleed flow to define the effects of tho bloods on dynatoic pressure character- 
istics. 
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4 . 1 , fl , 1 C’n vlty Geometry 


Profile and top view photographs of tho cascade for the 0,8H~rad (KWlog) 
Incidence, fl,l5H-m/n (approximately) Inlot voloolty series of tout point o nro 
shown In figures nnd in, Tho photographs nro arranged In ordor of Increasing 
cavitation number on onoh figure, nnd tho toot points, from tnblo 8, nro Indlontod 
undor onoh photograph. Mow In from loft to right, Counting hydrofollo from loft 
to right, tho profllo photographs nhow tho trnlllng edges of foil 8 nnd Colin 8, 4, 
find 0 nnd tho loading edge of foil 0, Tho lop view shows tho loading edge nuotlon 
luirfnoon of folh-i 8 (In part), 8, 4, R, nnd (1 (In part), Tho norlon of tout point n 
shown In tho photogrnphn 1 1 lustra tea tho nnturo of tho onnondo onvltntlon nnd In 
ropronontntlvo of othor norlon of polntn in that rnupoot. 

Tho dogroo of onvltntlon uniformity that oould bo ontnblbdiod o,i (ho various 
foils nnd spamvlsc on onoh foil onn bo noon In flguron 12 nnd 18, A 'air dogroo of 
uniformity was nooonipllnhod with tho nniiill cavities, but tho degree of nttnlnublo 
uniformity win not an good for tho lnrgor ouvlty size. Tho variation In onvltntlon 
pattornn In bollovod to bo Indlouttvo of vnrlntlonn In flow angle and/or voloolty. 

Tho pattornn woro mado fin uniform fin ponnlblo during tho sotting of onoh tost 
point through adjuntmont of tho oancndo nldo walls and dlschargo tall boards. 
Additional control Is possible for noncnvltntlng tost points through tho variation 
of timnol stop slot boundary layer removal flow (figure 7), but these slots cnvltntod 
when significant cavitation developed In tho cnscado nnd sufficient flow for two- 
dimensional cascade conditions was not believed to be established , The cavitation 
off the edges of the step slots can bo seen In figure 13, 

The observed cnscado nonuniformity nnd the occurrence of pressure surface 
cavitation at high Incidence ungles (0,23 rad, 13 deg) led to the conclusion that 
there were local variations in incidence angle and velocity within tho cascade, 
and that the average Incidence angle was lower than the Indicated geometric angle. 
It Is also probable that Incidence angle varied with cavitation number. It was not 
possible either to correct these conditions (through the design of a new cascade 
Inlet nozzle) or to measure the actual conditions (w'hlch would have required 
sophisticated Instrumentation such as a laser veloclmotcr) within the scope of the 
program. Cascade cavitation was generally unsteady, and this unsteadiness also 
varied from one blade to another. The unsteadiness seemed to be random and had 
no relation to cavity size. In light of these uncertainties, tho cascade cavity 
geometry tost data were considered to Indlcato only qualitative trends, and It was 
evaluated with the uncertainties In mind, 

A listing of tho test points, measured cavity lengths, nnd heights Is given 
In table 3, Cavity lengths for blades 2, 3, and 4- wete measured from the top 
view photographs at the foil mldspnn, where boundary layer effects would bo 
minimal. Blade 4 cavity length and cavity holght woio measured from the profile 
photographs. This blade was tho only one closo to a truo profllo In these photo- 
graphs, Cavity height was measured from tho blado suction surface In a direction 
perpendicular to tho pressure surface. All measurements were corrected for 
viewing angle. Cavity length measurements from the top view are missing from 
tho table In somo casos where tho cavity closure was under tho adjacent foil and 
could not bo measured. In such cases, the foil "4” profllo length measurement 
was used to establish length. 
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Tho boundary layor blood Riot vnlvoR wore closed to obflorvo what offoct 
boundary layor blood would have on tho cavitation pattern for throo points, Ha tod 
an to fit points No* J9 f 31 # and 41 In table 3* When the bleed vaJvos wero closed 
(no other notion wna taken), the euvltatlon pattorn changed significantly, The 
auction aurfaco cavltloa, which had boon from 8 to 2(1% chord In length, became 
very small misty patchoa, mid pressure aurfnoo onvttloa for mod. Tho proaauro 
surface cavltloa woro also very amall misty patches attest points No. 19 and OX* 
howovor, at tost point No. 43 thoy woro about 20% chord In length. Water velocity 
approaching tho cascado romnlnod rolatlvoly constant, whllo static prosHuro In- * 
creased, resulting In an Increase In cavitation numbor. Tho tncroaso In cavita- 
tion number explains the reduction In cavity size on tbj blade suction surfaces: 
howovor, the appearance of prossuro surface cavitation Indicates a change In flow 
angle (Incidence) or In uniformity as well. This shows a qualitative dependonco of 
the cavitation pattorn on blood flow. Cavitation was observed on tho slot Up at the 
lower cavitation numbers, but blade pressure surface cavitation did not occur In the 
Incidence range solected (0.28 to 0.36 rad, 16 deg to 20 deg) when the bleed valves 
were open. 

Measured cavity length and height data are shown plotted and correlated 
with predictions In paragraph 5.1.2. 


4. 1.6. 2 Cascade Instability 

A summary of the cascade leading edge transducer amplitude (No. 3 In flg- 
ure 6) and visual observations of cavity motion made during the tests are presented 
In table 4. An equipment malfunction that was unnoticed while the test points 
were being recorded caused the loss of tape recorded pressure data for points 
1 through 18. The leading edge transducer was the only one to show a significant 
amplitude during the tests. The upstream transducers (No. 2) and the trailing 
edge transducer (No. 4) were similar to transducer No. 1, typically showing a 
peak-to-peak amplitude of 0.1 N/cm 2 (0.15 psl) and never exceeding 0.2 N/cm 2 
(0.30 psl). The 0.1 N/cm 2 amplitude Is the level of normal tunnel vibration 
under noncavltatlng conditions. O-graph traces of transducers No. 2, 3, and 4 
at test points No. 28, 30, 36, 37, and 40 are shown In figure 14. These traces 
are typical and show both the highest and lowest amplitudes recorded at valid 
test points. It will be noted In table 4 that test point No. 19 has the highest 
“ffi? 1 ® 1 however * thls P 0 ^ and test points No. 31 and 41 wore conducted 
with the boundary layer bleed slot valves closed, which caused cavitation to 
occur on the blade pressure surfaces. These three "bleeds off" points are not 
typical and were not considered In the cascade Instability Investigation, 


In table 4 the visual observation of cavitation Is listed as unsteady or verv 
unsteady at almost all test points. Also, It will be noted that the visual observa- 
tlons do not agree well with the leading edge transducer amplitudes. The reason 

wLw 8 0 lS i the k} ade cavities almost never appeared steady In the sense of 
having a definite collapse point. The collapse point was continually ** ovlng In 

LnSl t? 0 ' 1 * Nor pos , 8lble t0 distinguish a smooth streamline 
separating the vapor cavity rrom the surrounding water. The cavity profile 

was generally Irregular and changed shape In a random fashion. To an observer, 
the Irregular shape and random motion gave the Impression of unsteady flow, 
even though the pressure oscillations wore no greater than for noncavltatlng flow. 
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Table 4. Cascade Dynamic Pressure Data and Observed Cavity Motion 



(1) Frequency Was Not Apparent on Oscillograph 









Amplitudes ond froquonclos recorded by the leading edge transducer aro 
plotted against cavitation number In figure ID, This figure rovonls no discernible 
trend of either amplitude or frequency. At cavitation numbers greater than 0,50, 
only normal vibration Is present, Hotwoon 0,05 nnd 0,30, both amplitude and 
frequency are random. Also shown In figure 15 for comparison Is the highest 
amplitude recorded by Acosta nnd Wade (6), Tho prosont loading odgo amplitudes 
aro of tho same order of magnitude as those measured by Acosta; howovor, Acosta 
recorded trailing edge amplitudes of the same order of magnitude as the lending 
odgo, while tho present trailing edge amplitudes novor changed significantly from 
normal vibration levels, as shown In the O-graph traces of figure 14, Tho roason 
may bo the proximity of the trailing odgo transducer to tho plenum tank, which 
effectively absorbed pressure oscillations. 
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Figure 15, Variation of Cascade Dynamic Pressure 
Data With Cavitation Number 
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Cavitation Is generally unsteady In the sense of continuous Irregular motion 
of the cavity, particularly In the collapse zone. However, unsteady cavitation does 
not necessarily mean that tho flow Is unstable. In experimentally determining 
whether a cavttatlng flow through a cascade or an Inducer Is stable or unstable 
the viewpoint adopted Is that the flow Is unstable only If upstream and downstream 
flowrates and pressures, or Internal flowrates and pressures, vary In some regular 
pattern. This definition Includes blade-to-blade propagation of cavitation and 
unisonous cavity oscillation as being unstable, and Is Intended to distinguish be- 
tween these two phenomena and cavitation noise generated by random Irregular 
cavity motion, which Is considered stable. 


Flow through the cascade was certainly unsteady, as witnessed by the 
Irregularities In cavity shape shown In the photographs and the visual observa- 
tions noted In table 4. However, from this evidence It Is difficult to characterize 
a given test point as stable or unstable In the sense of propagating cavitation or of 


I 


ssssii ?t-S”S 

edge pressure transducer tracos. 

In figure 15 the leading edge transducer amplitudes have been ^Mclodl^ 
three bands In an attempt to make a distinction botweon normal tunnel vibration, 
tXtton noise? and unstable Dow. Normal vibration levels wore easily Identified 
bv the amplitudes recorded at the higher cavitation numbers. The dividing point 
Kee^ and unstable flow was not easily Identified, and has been 

Strartir^ “d as the level separating the majority of the test point* .from 
the four test points showing the highest amplitudes. If any point were unstable, 
one would expect It to be one of the points with highest amplitude. A comparison 
of the visual observations, which are listed In table 4, with the four highest 
pressure amplitude points (No. 28, 36, 37, and 40) shows no Particular correla- 
tion between the dynamic pressure data and the observatlo ^‘ ^ st ^/^ t N 4 * 0 a 
Is listed as very unsteady, No. 36 as steady, No. 37 as pulsating, and No. 40 as 
steady. Figure 15, therefore, shows the only possible distinction between the 
experimentally stable and unstable test points. The points are discussed further 
and correlated with predictions In paragraph 5.2. 

4.2 INDUCER EXPERIMENT 

4.2.1 Facility 

TM 

The test facility used for the rotating Inducer Instability tests was P&WA 
FRDC test stand D-34, a schematic of which Is shown In figure 16. The 13-cm 
(5-ln.) Inside leg of the loop was used for all testing. 

The test rig is driven by a 260-hp do motor through a 2.33:1 ratio gearbox, 
providing output speeds to 9000 rpm. Rig speod Is automatically regulated to 
within 0.03% by control of tho dc output of the motor-generator set supplying 
drive motor power. The test loop uses demineralized water and has a heat 
exchanger for control of loop water temperature. A complete 
ts available for the deaeration of the test loop water. Loop pressures are con- 
trolled by either the Inlet or discharge accumulator system, depenulng on the 
Inlet pressure required. Tho Inlet accumulator Is usod for Inlet pressures 
above* 5 atmospheric, and the discharge accumulator Is used for lower than atmos- 
pherlc pressures. A schematic typical of both accumulator systems Is sjjovm ln 
figure 17. Inlot pressures down to approximately 2 N/cm2 ra.apsllarojosibic. 
Flowrate Is controlled by a motorized flow valve located In tho 13 cm (5 •) g 

and measured with a turbine-type flowmeter located downstream of the flow valve. 
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rho Inducor toot rig rlonlgn Ih shown In figure 1.8, The Inducer housing was 
nmdo of transparent acrylic and thuc permitted obRorvntlon of the Inducer cavita- 
tion, 'i'hc Inducer la overhung from Hr hearlnga In the Inlet line and Inlet flow Ir 
unobstructed. Figure 10 shows an overall view of the rig on the teat stand. 
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Figure 18, inducer Test Rig 
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Figure 10, Inducer Installed In 'lest Facility 
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4,2,2 inducers 


The two Inducers lifted In the tout program worn fabricated In tho Inducor 
load and stress project, NAHA Contract NAH,V112I0, an reported In (1). Tho 
inducers were identical except that the Ural Inducer had a radio) looking edge, 
while the aeeond Inducer wan a wept back at a cone angle of 0,28 rad (1(1 deg), in 
the manner ahown In figure 20. The conical cut wan moved axially forward the 
distance "Y" to provide for the fairing of the blade leading edge Into the hub. Thin 
Hweepback produced the lending edge contour ahown In figure 21, An the figure 
ahown, the awept loading edge Inducer had preaauro tap tubing Installed In tho 
bio don for the previous program. Tho blade aurfaee around the tubing- wan filled 
with epoxy to provide a am noth aurfaee, 

'J’he Inducer had three blades with Inlet tip blade angles of 0.14 rad (8 dog). 
The blading was helical to a solidity of approximately 1.0, at which point camber 
was applied exponentially, Tho swoop outback was within tho blade helical portion 
at all radii, so that blndo leading edge metal angle was the same for tho radial 
and swept loading odgo Inducers, Doth Inducers woro machined from titanium. 
Table 6 summarizes tho basic Inducer design parameters. 
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Figure 21. lncluccr With 0.28 rad (16 deg) I, ending FE 113407 

Edge Swecpback 


Table 5. inducer Design Parameters 


4,2,3 I nntru mentation 


Indmur porformonco Instrumentation was located no shown In figure 22. 
Inlnt total procure was measured approximately 80 om (30 In.) upntronm of (bo 
Inducer loading edge with n Kiel probn located nt tho Inlot plpo oontnrJlno, A 
tmvorflo wedge probe wan unod near tho 'idueor loading edge to monitor Inlot 
Uaw angle and flow velocity at a radlun »oor tho Inducer tip (8, (JO cm, 3,4J in.). 
I l f1u h f irgo ^otn * nnd static pressure nnd flow angle wore gathered with another 

£t V T? ^ f,R '° pr ^r* F r toflta, the probe waa net at a Slue 

that wan lepreaentatlvo of masa averaged totnl head, both Inlet and exit prohen 
were air calibrated before Inotnllntlon, and atatlo preoauroa readings were 
corrected accordingly. Static pressures at tho upntx*onm meaauromont station 
aud lm mediately upstream and downstream of tho Inducer wore measured with 
static wall taps. All pressures wore road from preolslon, a 0. 20% F.S. 
pressure gages. *’ 
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Figure 22, Inducer Pressure Instrumentation 
Locations 
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. D Tn am [ C p *? 8 l ur ? measurements wore taken with piezoelectric transducers* 
mountod flush with the Inside walls at the locations shown In figure 16. The ten 

iv«tn«, n8 H e r f e l°fu ec ! to pr , ovlde complote coverage of the Inducer-flow loop 

Output of the transducers was recorded on magnetic tape. Static pres- 

fnin+i * ^ UBh "m° un ted static wall taps were also gathered at each of the ten 
locations and read from gages. 


*Klstlor Model 601 H 
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A turblno-typo flowmotor located downstream of tho flow throttling valve 
In tho 13-om (G-ln.) log was usod to moasuro flowrate, Hoadout was on a digital 
counter, Wator temperature was measured with a chromol-nlumol thermocouple 
located upstream of tho Inducer Inlet, Inducer speed was moasurod by oloctronleally 
coimtlng tho output of a magnetic transducer that sensed tho passing frequency of 
a 00-tooth gonr on tho rig drlvo shaft. That signal was multiplied by a preset value 
and transmitted to a panel-mounted digital counter. 

Table 0 summarizes tho Instrumentation and the estimated maximum error 
associated with each rending, 

4,2,4 Tost Procedure 

Since the Inducers and the test rig that were usod In this program wero 
used In previous programs, the tests wero numbered consecutively from the 
previous programs. This program represented the eighth series of tests and 
each test Is designated 8. XX with the "XX" Indicating the test number for this 
program, 

4.2.4. 1 Performance Tests 

Prior to the start of Instability testing, noncavitattng performance tests 
were conducted on the radial leading edge Inducer to supplement the data reported 
In (1). Performance data were available for flow coefficients of 0.096, 0.090, 
and 0,084 and additional data were requ' red for 0.070. The available data had 
also been obtained with blade pressure Instrumentation on the Inducer, and the 
effects of the Instrumentation on performance were unknown. (The radial Inducer 
used In this program had no blade Instrumentation. ) 

The performance tests were conducted by setting a speed, flow, and Inlet 
pressure, and all other Instrumentation were recorded. Total pressure, static 
pressure, and flow angle were manually traversed at each of seven radial stations 
(spaced at equal area Increments) at the Inducer Inlet and exit. In total, three 
flow coefficients (0.096, 0.084, and 0.070) were explored. The speed was a 
constant 513.1 rad/s (4900 rpm) and Inlet pressure a constant 20.7 N/cm^ gage 
(30.0 pslg). 

4 . 2 . 4. 2 Accumulator E valuation 

Before the Instability testing was started, the effects of the test loop pres- 
sure control accumulators on the response of the total Inducer-test loop system 
was determined. Various combinations of "hard" and "soft" (wator or air In the 
bladder) accumulators were tried while pressuro oscillations were being generated 
with the radial leading edge Inducer. A schematic typical of the Inlet and discharge 
accumulator systems Is shown In figure 17. The Klstler transducers, mounted In 
the water loop at the locations shown In figure 16, were used to measure dynamic 
data during the time the different combinations of accumulators were being tried. 
Subsequent analysis of the dynamic pressure data showed it to be Independent of 
the accumulator configurations. Dynamic pressure amplitudes as high as 
±4 N/cm 2 (±6 psl) were generated at the Inducer discharge (location 10), 
diminishing to approximately ±1 N/cm 2 (±1.5 psl) at location 3. Negligible 
pressure amplitudes were measured at locations 4, 6, 6, and 7. The ampli- 
tude at location 8 was approximately ±1 N/cm 2 (±1.6 psl), and data at loca- 
tion 9 were obscured by high pressuro amplitudes at blade passing frequency. 
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Inducer Test Instrumentation (Continued) 
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a. 


Approximately 12 points wore selector! In this range to doflne the shape of tho 
head rlso vs cavitation number curve and to Identify tho range of unstablo opera" 
tlon. A camplote sot of stoady and dynamic pressure data wero taken at oach 
point, Tho cavitation was also obsorved and noted through tho acrylic housing 
with tho aid of a strobo light, 

4,2,5 Data Reduction 

Inducor performance data woro doflnod ns the head rise and offlcloncy 
botwoon tho upstream measurement station (Loop Station "8", PTUP In figure 22) 
and tho discharge measurement station (Loop Station "9", PTG2 In figure 22), 

The Klol probe at Station "8" was sufficiently far upstream that prerotatton did 
not offect Its readings, while the wedge probe at the Inducer Inlet (Station "9", 
PTG1 In figure 22) was effected by prerotation. Inlet conditions wore calculated 
from the data as follows: 


Upstream Total Pressure: 

Pg = PTUP (Corrected for Gage Height) 

Inlet Total Head: 

«1 - H 9 - Pg/p 

Inlet Absolute Velocity (Assumed Axial): V t = Q/Aj 

Inlet Static Head: 

h t = H t - Vi 2 /2g 

Inlet Tip Relative Velocity: 

V t '= '/Vf +U2 

Inlet Tip Relative Velocity Head: 

q t » • V‘2/2g 

Tip Cavitation Number: 

k = (hji - h v )/q* 

Inlet Tip Plow Coefficient: 

* - V u t 

Discharge conditions were calculated from the measured data at seven radial 
measurement stations, as foUows: 

Discharge Total Head: 

H d <= (corrected for gage height) 

Discharge Static Head: 

PSG2 1 

hjj » — - — (corrected for static pres- 
v sure and gage height) 

Discharge Flow Angle: 

0 d ■ BAB2J 

Discharge Velocity: 

v d = ^g(H d - h d ) 

Discharge Axial Velocity: 

V zd = v d 8ln 8 d 

Discharge Tangential Velocity: 

v ud “ v d 008 0d 
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Average ninehnrge Total Hoad: 


I 


Avoragv. Discharge Ideal Head: 


Average Head Coefficient: 


Avorage Efficiency: 


Y"t 


II. 
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(Assumes V u i = 0) 

Inlet total pressure, static pressure, and flow angle were measured at the inducer 
Inlet tip and were reduced as follows: 


Inlet Total Head: 

PTP1 T 

Hi = — p — — (corrected for gage height) 

Inlet Static Head: 

PSG1J 

hi = — jg — (corrected for gage height) 

Inlet Flow Angle: 

= BAB1J 

Inlet Velocity: 

V l= :l/2g (H t -h t ) 

Inlet Axial Velocity: 

V zl ° V l sln 01 

Inlet Tangential Velocity: 

V ul = V l cos 0 l 


Nomenclature Is defined In the Appendix. Measurements and calculations were 
performed using the English system of units and were converted to the SI for 
presentation In this report. 

The magnetic-tape-recorded, piezoelectric, transducer data were reduced 
with the aid of a spectrum analyzer. A block diagram of the system used Is 
shown In figure 23. The Input (test data on the magnetic tape) was Input to 
an amplifier and multiplied by some factor N. This amplification was necessary 
to ensure that the data entering the spectrum analyzer were within the analyzer 
Input level limits. The spectrum analyzer contained an internal 50 Hz low-pass 
filter with a 0. 16 Hz noise bandwidth, A 10-sec record length of transducer data 
was put Into the spectrum analyzer’s memory after passing through this filter. 
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Thin Information wan converted to n frequency vn amplitude output with a resolu- 
tion of 0,1 IIz, After panning through an atlonuator to null out the offoots of tho 
first amplifier, the output was scaled for plotting purposes and thon plotted on 
paper by a Moseley X-Y recorder, Tho resultant plots had a frequency scale of 
0 to 00 II z, and are accurate at all frequencies nbovo 0,3 Hz, which Is the lower 
limit of the data response and reduction system. 

Input to the amplifier and output of the scaler wero monitored with oscillo- 
scopes, as shown In figure 23, to onsuro corroct system operation. Prior to the 
tapod data being Input, a calibration signal of known amplitude and frequency was 
Injected Into tho system and proper functioning of tho entire setup was chocked. 
The tapo reproducer was Identical to the recorder used to gather the data at tho 
test stand, thus ensuring that what canto off the tapo was Identical to what wont on. 



Figure 23. Data Tapo Analyzer System FD 63040 

4.2,6 Test Results 

A list of Inducer tests and test point conditions Is given In table 7. The 
first three tests (8.01, 8.02, and 8.03), were noncavltatlng performance tests 
for tho radial leading edge Inducer; tests 8. 04 through 8. 08 were radial Inducer 
stability tests; anti tests 8.09 through 8. 11 were swept Inducer stability tests. 
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Tablo 7. Summary of Inducor Tost Points 


tfpocd 

rnd/s rpm 


Hadlal i .ending I'Mgo Inducor 



Inlet Total 

Water 

Plow 

Pressure 

Tomporaturo 

cm 3 /a gpm 

N/cm 2 pslg 

«K H F 


01 

513.1 

4900 

02 

613.1 

4900 

03 

513.1 

4900 

,04 

513.1 

4900 


513.1 

4900 


513.1 

4900 


613.1 

4900 


513.1 

4900 


513.1 

4900 


513.1 

4900 


613.1 

4900 


513.1 

4900 


513.1 

4900 


513.1 

4900 


613.1 

4900 


513.1 

513.1 

513.1 

613.1 

513.1 
513.1 
513.1 
513.1 
513.1 
513.1 

513.1 

513.1 

613.1 

613.1 

513.1 
513.1 

513.1 

613.1 

613.1 

513.1 

613.1 
613.1 


4900 

4900 

4900 

4900 

4900 

4900 

4900 

4900 

4900 

4900 

4900 

4900 

4900 

4900 

4900 

4900 

4900 

4900 

4900 

4900 

4900 

4900 


85,940 

85.940 

86.940 

86.940 

85.940 
85,940 

85.940 

86.940 

85.940 
85,940 
85,940 
85,940 

80,260 

80,260 

80,260 

80,260 

80,260 

80,260 

80,260 

80,260 

80,260 

80,260 

66,890 

66,890 

66,890 

66,890 

66,890 

66,890 

66,890 

66,890 

66,890 

66,890 

66,890 

66,890 


1362 

1362 

1362 

1362 

1362 

1362 

1362 

1362 

1362 

1362 

1362 

1362 

1272 

1272 

1272 

1272 

1272 

1272 

1272 

1272 

1272 

1272 

1060 

1060 

1060 

1060 

1060 

1060 

1060 

1060 

1060 

1060 

1060 

1060 


14.6 

11.2 

8.0 

3.8 

1.0 

0.3 

-0.3 

- 1.0 

- 2.1 

-2.4 

-2.5 

-4.9 


( 'ommontH 


80.0 302.4 85 Noncavltntlng Por- 

formnneo Test 

30.0 304.1 88 Noncnvltntlng Per- 

formance Test 

30.0 306.7 91 Noncavltntlng Per- 

formance Tost 


21.2 

16.2 

11.6 

5.5 

1.4 

0.4 

-0.46 

-1.4 

-3.0 

-3.5 

-3.7 

-7.1 


302.4 

303.0 

304.1 
304.7 

305.2 
305.2 
305.2 
305.2 
305.2 
305.2 
305.2 
305.2 

302.4 

302.4 

303.6 

304.7 


-V/# V 
-2.0 

x ♦ u V vv f 

-2.9 306.2 

90 

-2.4 

-3.6 306.2 

90 

-3.1 

-4.6 306.2 

90 

-4.8 

-6.9 305.2 

90 

-6.0 

-8.7 306.2 

90 

14.8 

21,6 300.6 

82 

8.1 

11.7 303.0 

86 

4.7 

6,8 303.6 

87 

1.2 

1,7 304.1 

88 

0.8 

1.2 304.1 

88 

-0.6 

-0.9 304.7 

89 

-1.8 

-2.6 304.7 

89 

-2.4 

—3.6 306.2 

90 

-4.4 

-6.4 306.2 

90 

-5.7 

-8,3 306.2 

90 

-6.9 

-8.6 305.2 

90 

-6.1 

-8.9 305.2 

90 


Test loop water con- 
tained excessive 
amount of air 
(0.0148 m 3 , 4 gal, 
registered during 
deaeration test, as 
compared with 
0.0037 m 3 , 1 gal, 
for all subsequent 
runs). 


Discharge accumu- 
lator system de- 
veloped small line 
leak during test, 
thereby causing dif- 
ficulty In holding- 
low pressure test 
points steady. 


Normal 
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Table 7. Summary of Inducer Teat Points (Continued) 


Speed 

Flow 

n , 

Inlet Total 
Procure 

Water 

Tomporature 

rad/s rpm 

cnr/fl gpm 

N/em* palp 

°K °P 


Comments 


Hadlal Leading Edge Induaor 

8.07 

513. 1 

4900 

85,940 


1 

15.0 

21.7 

302.4 

85 



613,1 

4900 

85,940 

136 

2 

7.9 

11.6 

304.1 

88 



513.1 

4900 

85,940 

136 

2 

4,5 

6.5 

304.1 

88 



613.1 

4900 

86,940 

1362 

1,0 

1.5 

304.1 

88 



513.1 

4900 

85,940 

1362 

-0.5 

-0.7 

304.1 

88 

Repeat of tost 8.04 


613.1 

4900 

85,940 

1362 

-1.1 

-1.6 

304.1 

88 

with water deaerated 


613.1 

4900 

85,940 

1362 

-1.8 

-2.6 

304.1 

88 

to stand limit. 


613.1 

4900 

85,940 

1362 

-2.4 

-3.6 

304.1 

88 



613.1 

4900 

86,940 

1362 

-3.1 

~4#5 

304.1 

88 



513.1 

4900 

85,940 

1362 

-4.8 

-7.0 

304.1 

88 



513.1 

4900 

85,940 

1362 

-5.6 

-8.2 

304.1 

88 


8.08 

513.1 

4900 

80,260 

1272 

14,9 

21.6 

304.1 

88 



513.1 

4900 

80,260 

1272 

1.1 

1.6 

304.1 

88 

Repeat of test 8, 05 


513.1 

4900 

80,260 

1272 

-2.6 

-3.8 

304.1 

88 

at lower pressures 


613.1 

4900 

80,260 

1272 

-4.0 

-6.8 

304.1 

88 

with properly func- 


613.1 

4900 

80,260 

1272 

-5.6 

-8.2 

304.1 

88 

tioning discharge 


513.1 

4900 

80,260 

1272 

-5.8 

-8.4 

304.1 

88 

accumulator system. 


513.1 

4900 

80,260 

1272 

-5.8 

-8.4 

304.1 

88 


0.28 rad (16 deg) Swept Leading Edge Inducer 


15.0 

21,7 

8.1 

11.7 

4.6 

6.7 

0.3 

0.4 

-1.7 

-2.4 

-2.2 

-3.2 

-2.5 

-3.7 

-4.3 

-6.3 

-5.0 

-7.2 

-6.2 

-9.0 

14,8 

21.5 

8.1 

8.1 

4.7 

6.8 

1.2 

1.8 

0.3 

0.4 

-1.7 

-2.5 

-2.6 

-3.6 

-4.1 

-6.0 

-4.8 

-7.0 

-5.9 

-8.6 

-7.2 

-10.5 












Tablo 7. .Summary of inducer To fit Points (Continued) 


Tost 

No, 

Spood 

rad/s rpm 

Flow 

cm 3 /a gpm 

Inlet Total 
Pressure 
N/em 3 pslg 

Water 

Temperature 
°K °F 

Comments 

8,11 

013.1 4000 

513.1 4900 

013. 1 4900 

013.1 4900 

513.1 4900 

(10,890 1000 
00,890 1000 
00,890 1000 
00,890 1000 
00,890 1000 

14.8 21.0 

8.0 11.0 

4.0 0.7 

1.3 1.8 

0 0 

299.7 80 
301.2 83 

301.8 84 
302.4 85 
303.0 80 

Induce** housing split 
midway through test, 
forcing termination. 


Indue or head rise coefficient and efficiency data are plotted against How 
coeflldent In figure 24. The radial Inducer data are from tests 8.01 through 8. 03 
and the swept Inducer data arc from the previous program reported In ( 1 ). Radial 
Inducer tost data, with and without blade Instrumentation, were found to agree 
almost exactly with provlous data at the two flow coefficients that were re- 
peated (0*096 and 0.084). It can therefore be concluded that the pressure 
Instrumentation that was on the radial Inducer blades for the previous tests of (1) 
and also on the swept Inducer tested In this program, had a negligible effect on ’ 
overall performance. 




Figure 24. Noncavltatlng Inducer Performance DF 91153 
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Tents B, 04 through 8,08 woro rncllnl loading odgo Induuor ntablllty tostp., 

Thn Indueor wan operated at a rotatlvo apood of 018 rad/a (4000 rpm) and flown 
of 80,940, 80,860, and 60,800 omVfl (I860, 1878, and 1060 gpm), Tnlot tip 
flow eoofflclontn at thono conditions woro 0,090, 0,084, and 0,070, respectively, 

In all canon thn objective wan to drop Inlot pronnuro to tin* point whom hond 
rlno foil off by at loant 80%, Originally throo tontn woro to bo conductod, ono at 
oaoh flow coofflolont# but tho doaoratlon attempt prior to toot 8. 04 represented 
tho flrnt tlmo that tho components of tho doaoratlon nyntom had boon operator! for 
Romo tlmo, and novoral nmall problomn showed up In tho prooorluro, An a result, 
tho loop wator oontalnod more dissolved air than that which wan present during 
tho later tontn, (Refer to the comments on table 7. ) In addition, during tent 
8.04 a loose coupling wan noticed In tho drive train after a cavitation number of 
0, OSS wan reached, and the tost wan terminated nbovo tho hond falloff point. For 
those reasons, tent 8. 07 ropoats the conditions of tost 8. 04. Tost 8. 08 was 
conducted to supplement tost 8. 05 at tho low ond of the cavitation number scale 
because the discharge accumulator system developed a small leak on tho surge 
tank side during tho first test at this flow coefficient. This leak mode It difficult 
to hold precise test points at low Inlet pressures. 

Cnvltatlng performance data for the radial leading edge inducer are shown In 
figure 25. Regions of observed test loop surging . re Indicated on the curves. It 
was Impossible to reach head falloff at the 0. 090 flow coefficient as tho Inducer 
loop system became uncontrollable at a cavitation number of approximately 0. 030. 

Large flow fluctuations of approximately 38,000 cm 3 /s (600 gpm) at approximately 
two cycles per minute prevented operation nt lower cavitation numbers. The 
results of tost 8.04 have been omitted from those curves since the loop water 
contained four times the amount of air that was present In all subsequent tests, 
and this excess air casts doubtB on the usefulness of the data for test to test 
comparisons. Surge regions were characterized by a cyclic chugging sound and 
by alternate lengthening and shortening of the tip clearance cavitation. Blade 
surface cavitation could not be clearly seen through the tip clearance cavitation. 

Two regions of surging were noted in the radial leading edge tests. The first 
region started at a cavitation number of approximately 0. 07 and completely 
ceased at a lower cavitation number before the start of the second surging region. 

The first region was always relatively mild as compared with the second. The 
second surging region ceased just prior to head falloff for the 0. 070 and 0. 084 
flow coefficients; however, as stated previously, It became very severe at the 
0. 090 flow coefficient and was accompanied by large test loop flow fluctuations. 

The 0. 28 rad (16 deg) leading edge Inducer was tested In runs 8. 09 through 
8. 11 In the same manner as the radial inducer. These test points are also summarized 
in table 7. Again it was impossible to reach a head falloff condition at 0. 090 flow 
coefficient. The same fluctuations that occurred with the radial inducer occurred 
again below a cavitation number of 0. 027. During the final test, 8. 11, the acrylic 
Inducer housing split and forced the termination of the run after the fifth data point.. 

The housing had been experiencing very heavy cavitation damage near the blade 
leading edge area during the course of the program. Cavltat. j performance data 
for the swept leading edge Inducer are shown In figure 26. Regions of observed 
test loop surging are again Indicated on the curves. The swept inducer exhibited 
only one surging region that occurred over approximately the same range of 
cavitation numbers as the second region of the radial leading edge Inducer. Surging 
with tho swept Inducer was milder than that which was observed with the radial 
Inducer. 
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Figure 25, Cavltatlng Performance, Radial Loading DF 91164 
Edge Inducer 
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Figure 26. Cavltatlng Performance, 0.28 rad DF 91165 

(16 deg) Swept Leading Edge Inducer 
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Figure 27, Normalized Ca vita ting Performance, 
Radial Leading Edge Inducer 
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Figure 28. Normalized Cavltatlng Performance, 
0, 28 rad (16 deg) Swept Loading Edge 
Inducer 
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An oscillograph trace of tho tape recorded dynamic pressure data for font 
B, 04 at Itn lowest pressure point (k - 0.008) in shown In figure 20. Thin trace In 
typkm! of tho relatively high nmpllturlo measured inntohllltlon found in nil tontn, 
nnd It won unod to gain nn appreciation of tho form of tho pronnuro signals and 
of tholr rolntlonnhlp to ono nnothor. Tho transducer ponltlonn nro shown on 
I'lguro it) Tho discharge transduce™ (10 nnd 1) nhow tho highest amplitudes nnd 
n Ini 1 1 n r oxc*opt tlint No* 1 0 1 which In nonvnr the Indueor, nhown n nl| r nl.flr*nnf 
H2 Hz onmponont (rotational npood) In spite of tho 00 Hz low pnnn tutor. No. l 
hnn n olonnor signal (with regard to tho donlrod cavitation Instability) nnd wnn 
selected nn tho representative trsnaduccr for dntfi correlating purposes, No. 8, 
In tho Inducer Inlet lino, nhown n measurable nlgnnl nt tho nnmo low frequency 
mid In phnno with tho signal from No. 1 nnd nine him nn 82 Hz component. Thin 
nlgnnl wnn nine considered to ndor|untoly roprenont tho Inntnblllty. TrnnHducor 
No. 0| nt tho Inducer Inlet, hud Its dntu completely obscured by blndo nnsHlnu 
frocjuonoy (24(1 Ilz). H 


Spectrum analysis plots of tho discharge transducer signal (No, 10) over a 
0. 0 to on Ilz) range nro shown In figures DO through ,17. Ono figure shows the 
p.lotfl for nil tho teat points of a given tost* Tho plots gonernlly show that the 
dynamic presnure signal Is made up of a woll-doflned, fundamental frequency with 
harm on? eo. The harmonics frequently are of higher amplitude than tho fundamental 
frequency. Note that amplitude scales are not tho same for all of the plots. 


, transducer (No. 8) spectrum plots were uniformly similar to thoso 

from the discharge transducer except that they wore of lower amplitude. Several 
of the inlet plots that had significant amplitudes are shown In figure 38. 


Since the discharge transducer signal was found to adequately represent the 
system dynamic pressure data, the amplitude and frequency of Its fundamental 
component were plotted against cavitation number for each flow coefficient. These 
data are shown for the radial leading edge Inducer in figure 39 and for the swept 
leading edge Inducer In figure 40. Amplitude data are plotted even where the 
spectrum Is apparently composed primarily of noise, whilo frequency data are 
plotted only where a significant frequency could be Identified on the spectrum plot. 
For points where the pressure amplitudes were low, the identification of a 
fundamental frequency requtrod engineering Judgment) figures 39 and 40 should 

b t. e I^ Uat , Gd Wlth ^ ls mlnd » T h e frequency curves are quite smooth, howover, 
which lends considerable confidence to the accuracy of the spectrum data Inter- 
pretation. The location of the observed surge regions Is indicated on the figures. 
Surge regions wore characterized by visual and audible flow oscillations, as de- 
scribed previously. 

Both tho radial and tho swept Inducers have Identifiable discrete pressure 
oscillations at relatively high cavitation numbers (0. 16 to 0. 23). The amplitudes 
tend to gradually Increase as cavitation number Is lowered until the seeond observed 
surge region of tho radial tnduccr and tho only surge region of the swept Inducer 
(figures 26 and 26) ts reached. At this point there Is a sudden Increase In pressure 
amplitude. Thorc Is no particular correlation between the pressuro amplitude data 
and the first observed surge region oi tho radial Inducer; however, these surges 
were relatively mild and often difficult to detect so pressure amplitudes would not 
be expected to be particularly high. Pressure amplitudes were generally con- 
siderably lower for the swept inducer than for the radial. 




Figure 29. Oscillograph Tracing Typical of Inducer FD 63039 
Instability? Test 8,04, 513,1 rad/s 
(4900 rpm), 5= 0.090, k = 0.038 
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Figure 30. Inducer Discharge Dynamic Pressure Spectrum Plots; Radial Leading Edge, FD 62545 

0 = 0.070 (Test 8.06) 
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Figure 40. Variation of Inducer Discharge Dynamic 
Pressure Data With Cavitation Number. 
0.28 rad (16 deg) Swept Leading Edge 
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Frequencies tend to deerenne gradually to a captation number of approxi- 
mately O.Of) for tho radial Inducer, and than to dnerenne more rapidly, Tho flrat 
aurge region occurred when I'roquonolen worn In tho H to Jo Ilx range, Froquenclea 
wore generally higher for tho swept Induoor and a rapid decrease with cavitation 
number atnrta at a lower cavitation number (0, (in to o, oil). 

Considering tho variation of tho data with flow coefficients, tho radial 
Induoor ahown a definite amplitude Ineronao an How ooofflclont docroaaon (at 
cavitation numbora where tho measured amplitudes wore reaaonably largo); tho 
awopt inducer data are limited and do not ahow a definite relationship. The 
radial Inducer aliowa no particular relatlonahlp of frequency to How coefficient' 
while the awopt Inducer aliowa a definite reduction of frequency with flow 
coefficient, 

A sight glasa htul boon Installed In tho tost loop at tho Inducer Inlet accumu- 
lator (No. 7 In figure 16) for tho accumulator evaluation toata. The lovol of water 
In this glass was noted for tests 8.00, 8.08, and 8,09 and a froo surface was ob- 
served to suddenly form at an Inlet pressure of approximately -4.5 N/cm 2 (-6.5 pslg, 
Si 0.06). The formation of this surface was attributed to the flashing of air 
from solution and Is therefore Indicative of the achieved degree of water deaeration. 
Since the deareatlon procedure was the same for all tests, It can be assumed that 
the free surface generally formed at the same Inlet pressure for all tests. This 
Inlet pressure approximately corresponds to the cavitation number at which the 
second surge region for tho radial Inducer and tho only surge region for the swept 
Inducer were noted. Hence, there was a sudden change In the test loop system 
that coincided with the observed surge regions. 

The general effects of air In the loop system on the pressure signals can be 
evaluated by considering the spectrum plots for test 8, 04 (figure 33) with those 
for test 8.07 (figure 34), Both tests were run with the radial Inducer at a flow 
coefficient of 0. 090, but test 8. 04 was not considered satisfactory because tho 
water Initially had approximately four times as much air In solution as did test 
8, 07 and all other tests. Spectrum plots for test 8. 04 show noticeably less well 
defined signal frequency components than the plots for test 8. 07, and test 8.04 
had relatively high amplitude, low frequency oscillations at a cavitation number 
of 0,038, whereas test 8.07 did not at a nearly comparable cavitation number of 
0,039, The amplitude and frequency of the signals are plotted against cavitation 
number In figure 41. The data for the two tests are similar except that test 8.04 
tends to show slightly lower frequency oscillations and Its pressure amplitude 
Increases at a higher cavitation number. 

Inducer Inlet fluid tangential velocities, which were measured near the 
Inducer tip, are shown plotted against cavitation number In figure 42 for the 
radial leading edge Inducor. The probe was located at a radius of 8. 66 cm (3.41 In. ) 
and was 0, 66 cm (0 26 In. ) upstream of the leading odge. Tip radius wub 8. 90 cm 
(3, 50 In. ). Tangential velocity generally decreased with decreasing cavitation 
number and reached zero at a cavitation number of 0. 06 to 0. 07 for all flow 
coefficients, the same range of cavitation numbers at which the first surge region 
was observed. Although not plotted, fluid through-flow velocity decreased along 
with tangential velocity. Tho rate of decrease with cavitation number was 
gradual to a cavitation number of approximately 0. 10, at which point velocity 
decreased more rapidly. 


73 


Ml 


Ml 


J 1 «.0 


3,U 


£ 


1.0 


0 


I* 11,0*10 


Heavy i Q 

Surging [ \ 

< ■'» V 


W O’ 




0 0.U2 0.0* 0,00 0.08 0.*0 


I i"H 

i i IMI'I (Ni'iMliMi'raitdi 
♦ '♦Ml/ 1 1 lead, Hall 


"Q <>' ^ (jQ> 

0. iT "ala 0 16 ♦ 0 22 M4 


CAVITATION NUMU1.H. k 



Figure 41. Comparison of Relatively Deaerated With DF 91160 
Nondeaerated Dynamic Pressure Data, 

Radial Leading Edge 
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Figure 42. Variation of Inlet Tip Prerotation With DF 91161 

Cavitation Number, Radial Leading Edge 
Inducer 


Tangential velocity data for the swept inducer are shown In figure 43. 

The probe was at the same radius but the sweepback of the blades resulted in 
the probe being 2. 20 cm (0. 87 In. ) upstream of the leading edge tip. Data were 
not taken for all the swept Inducer test points because of the time required and 
the observed high rate of cavitation damage to the Inducer housing. The data 
show the same trend as the radial data, decreasing to zero at a cavitation 
number of 0. 06 to 0. 07. Velocities were lower for the swept: Inducer than for 
the radial, but this may be at least partly explained by the probe location, which 
was approximately three times further upstream of the leading edge. 

Test results are discussed further and correlated with predictions i.n 
paragraph 5. 3. 
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Figure 43. 


v a riatton of Inlet Tip Prerotation Wit 
Cavitation Number, 0.28 rad (16 dee) 
Swept Leading Edge Inducer 8 
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SECTION 5 

ANALYSIS AND DATA CORRELATION 


Three sots of experimental data woro correlated with predictions; moasurod 
hydrofoil cascado suction surfaco cavity goomotry data, moasurod cascade tunnel 
pressure oscillation data, and measured Inducer loop Instability data. The cascade 
geometry data woro correlated against predictions made through the uso of an In- 
ducer Internal flow analysis, cascade tunnel dynamic pressure data were corre- 
lated with tunnel system model predictions, and Inducer loop dynamic pressure 
data were correlated with Inducer system model predictions. The analyses and 
results of each of the correlation efforts are reported In the following subsections. 

5. 1 CAVITATION CAVITY GEOMETRY 

5.1.1 Analytic Model 

The analytic cavity model was formulated In the Inducer flow analysts re- 
ported In (1) and Is based on the concept of a distinct vapor cavity on the blade 
suction surface which displaces an otherwise incompressible liquid. The liquid 
and vapor are assumed to be separated by centrifugal effects both In the meridional 
and blade-to-blade directions Such a cavity Is Illustrated In figure 44. The co- 
ordinates of the vapor-ltquld Interface (which defines the cavity) are calculated 
to the point of maximum cavity height In an Iterative process where the flow forces 
associated with the cavity are calculated and equated to the change In momentum 
parallel to the blades, and the radial streamlines are located to satisfy radial 
equilibrium. The cavity Is assumed to collapse linearly from the point of maxi- 
mum height In a manner that Is based on empirical data. The profile of such a 
cavity Is shown In figure 46. Significant dimensions are: 

c = Length from blade leading edge to point of maximum 
height along the blade chord 

b = Length of collapse region along the blade chord 

h = Cavity height normal to the blade chord measured from 
the suction surface 

1 = Length from blade leading edge to end of cavity along 
the blade chord 

i w = Angle of the cavity collapse (or wake) relative to the 
blade chord. 

The cavity model Is similar to the models of Stripling and Acosta (19) and 
Jalcobsen (24) who’ applied conformal mapping to obtain an exact solution for two- 
dimensional, Ideal flow. These models do net account for real effects, however, 
which occur In viscous, three-dimensional flow and which are accounted for In 
the cavity model of (1). Viscous effocts Influence the cavity through boundary 
layer blockage and drag forces. Three-dimensional effects consist of (1) centrif- 
ugal forces, (2) streamline relocation, and (3) changes In pas sag' width and height. 
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The analytic formulation of tho cavity model Ib given In (1) and It Ib In- 
cluded In the Inducor hydrodynamic computer program llstod In Volumo II of 
that roport. Correlations of tho model's predictions with test data and refine- 
ments to the model that wore made In this program are reported In the following 
subsections. 

5.1.2 Correlation of Cavity Model Predictions and Refinements of Model 

Cavity model predictions were correlated with predictions from the equa- 
tions of Stripling and Acosta (19), with cascade test data that were reported In 
Section 4. 1, and with Inducer test data from (1) to define areas of the model 
that required refinement. In this correlation effort, the dimensions shown on 
the cavity profile of figure 45 were assumed to define the cavity. Since Stripling 
provides an exact solution of the cavity profile to the point of maximum height 
for Ideal two dimensional flow with zero blade thickness, their predictions were 
correlated with cavity model predictions of height (h) and length to maximum 
height (c) for the same zero blade thickness, Ideal flow cascade. Cavity model 
predictions for such a cascade will differ from those for a real flow and real 
cascade but the degree of correlation provides an Indication of the cavity model's 
general accuracy. 

The data correlations and model refinements are discussed In the following 
paragraphs for each of the significant dimensions. Cavity model predictions of 
height were found to be reasonably accurate, length to maximum height was In 
error but correctable through adjustment of the empirical collapse function, and 
overall length predictions correlated reasonably well with the available test data. 

5. 1.2.1 Cavity Height 

Comparisons of cavity height predictions (dimension h In figure 45) with 
predictions from Stripling and Acosta are shown In figure 46 for the hydrofoil 
cascade and In figure 47 for an Inducer cascade. Actual height Is divided by 
blade spacing [7.6 cm, (3.0 In.) for the hydrofoil cascade and 18.6 cm (7.3 In.) 
for the inducer cascade] and plotted against cavitation number. The Inducer 
cascade Is geometrically similar to the tip streamline geometry of the Inducer 
that was tested as reported In (1). Cavity model predictions for both cascades 
were generated for an Invlscld, two-dimensional flow and zero thickness blades 
to provide a direct comparison with exact Stripling and Acosta predictions. 

Since the Stripling and Acosta model is an exact solution, the degree of correla- 
tion Is an Indication of the model's accuracy for the selected cascades and flow 
conditions. Figure 46 shows the cavity height correlation for the Ideal hydrofoil 
cascade at a single Incidence angle (0.28 rad, 16 deg). The correlation between 
the two prediction models Is good over the entire range of cavitation numbers. 
Figure 47 shows the height correlation for the Ideal Inducer cascade for three 
Incidence angles. At high cavitation numbers (1.0-0. 2), the correlation for all 
Incidence angles Is good and at low cavitation numbers the cavity model predicts 
cavity heights approximately 10% lower than those from Stripling. The low pre- 
diction la consistent for all Incidence angles. The Indicated degree of correlation 
for both cascades Is considered satisfactory and the wide difference between the 
blade stagger angles and Incidence angles for the two cascades Indicates that the 
cavity model Is applicable to a wide range of geometries and flow conditions. 
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Figure 46. Comparison of Cavity Height Predictions DF 91076 
With Exact Predictions; Ideal Hydrofoil 
Cascade l = 0.78 rad (16 deg) 
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Figure 47. Comparison of Cavity Height Predictions DF 91092 
With Exact Predictions, Ideal Inducer 
Cascade 
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Cavity model predictions of cavity height are compared with the hydrofoil 
cascade tent data In figure 48, 'i'ho measured and predicted cavity height data 
In this figure are from the auction surface of the hydrofoil, Tost data scatter la 
attributed to difficulties) oxporleneod In taking accurate moaauromonta from the 
photographa , to the fact that cavitation svna unatoady for noarly all teat points, 
to observed blado-to-blado nonunlformlty, and to possible variations In Incidence 
angle with cavitation number (duo to choking of the boundary layer romoval slots). 
Measured data were at Indicated Incldonco angles of 0,28, 0.81, a;.d 0.36 rad 
(10, 18, 20 deg) and the prodlctod data are shown for 0,28 and 0,17 rad (10 and 
10 deg). Cavity height predictions fall within tho band of the measured data 
scatter In spite of tho lower Incldonco angles for tho predictions. This result Is 
consistent with observed pressure surface cavitation at relatively high Incldonco 
angles and the conclusion that actual cascade Incidence angles were loss than tho 
value measured from tho geometric settings. 


(ncldcnct Velocity 



Figure 48. Comparison of Cavity Height Predictions DF 89778 
With Cascade Test Data 

The good correlation between height predictions from the cavity model and 
the oxact predictions from Stripling and Acosta for Ideal two-dlmons tonal flow 
led to the conclusion that tho model's predictions of cavity height are satisfactory. 
The hydrofoil cascade test data were not sufficiently accurate for quantitative 
correlation and refinement of tho model (as was originally Intended) but the 
model's predictions are In qualitative agreement with the data. This agreement 
also tends to substantiate the accuracy of the holght predictions, 

6. 1,2.2 Cavity Length to Maximum Height 

Length to maximum cavity height (dimension c In figure 46) Is precisely 
predictable, by the method of Stripling and Acosta, for an Ideal two-dimensional 
flow. Thus, while this length was not measurable In either the hydrofoil cascade 
or the Inducer test program of 0, a comparison of Idoal flow predictions from 
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m^STcouCr. 101 '’ m MMiUmn> ' n, " mtt0a 0 ‘ th0 '"" uo “ r “ovlty 

upn f ] , 0n, ?M t0 mnx ! nu, m holfSht predictions nro shown In fig- 

m o 40 loi the Irlonl hydrofoil cuh carlo nntl In figure 00 for tho Idonl tnrluoor cascade 
Jho bond that Is Indlcniorl for tiio Inducor program cavity model predictions is n * 
result of the ilnlto difference solution and It doflnos the limits of the length nretlle 
tlon at a Riven cavitation number. The predicted band wa^not 4own ln the height 
in^ f +i C ^ f nS ^ ur( r w 47 48 because height changes more gradually than length 

and die bond Is small. Tho cavity model's prediction of lontfh to maxfmum hoSt 
i the Ideal hydrofoil cascade generally shows tho samo trend with cavitation 
number as the Stripling and Acosta prediction but tho model's prediction Is longer 
“ '«* s “tP*l»S ™> d Acosta solution* Cnvtty modor p rotUc®o„s?orth° 

Ideal Inducer cascade are considerably longer than tho Stripling predictions at 
high cavitation numbers but they approach the Stripling predictions as cavitation 
7* th f *T T oquai at 8 numbor of^pproxlmatoly 

m-cdlot^mSi='.i 0r ?! °f v ! ty leaeth wltl > Incidence Is similar for both 
prediction models, These rosults Indicate that the cavity model generally ore- 

fmnnJfntt? e f 2T actU ,^i° ngths to maximum height. This Is considered to be an 
L def c en w. n the ™ avlty model and refinements were considered to 
improve the correlation. Tho calculations of cavity height and length to maximum 
height arc integral with the inducer internal flow analysis. hSwevef ^d^ S 
cments in the calculations would require major rework of that analysis A 
s mpler, more easily implemented, approach was taken wherein the cavitv 
collapse length (b in figure 45) was adjusted to make predictions of overall lemrfii 

fimofh ™ lth th< T test ^ ata * Thls approach was possible because the trend of the ** 
length to maximum height predlct'ons was correct and the cavity model's predic- 
n of length to maximum height was not enough longer than fhp r. nr i 

Acocta exact length to require collapse angles^* 

1.57 rad (90 deg). The error In cavity volume which Is thZ^ .m?” 

?h P nXXVubsimtonr nta|ie ° f tht ’ °' ,er "“ V ° 1Ume aa wU1 *» “Pirated In 

5. 1. 2. 3 Cavity Collapse 

it Th ?, cavlty c °Uapse model that was reported in (]) assumed that th P nnuit,, 

The empirical collapse angle function was reexamined in this orogr-im In n» 

cS °enS’ r ThXork r f ta°‘‘ T tW ? P rc,llcU " 1 «“• measured inducor overall 
cavity length. The work or other investigators was reviewed In r1r>w™i»„ , r « 

assumed form of the collapse luncion was cor,*ct AblSIt and Kl noT^stn L 

XsSv s^rnUnr ta n‘‘ 0n | b ^ l " d backwa ‘’ d step, a situation ttal ?s 

pnysicaiiy similar to flow in the wake region of a cavitation nnvttv ci u 

the experimental configuration that was used by Abbott aid Kline along with 
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to bo linear (n noooRRnry approximation In Wow of the lack of nny dntn concerning 
tho notunl shape) from tho Rtop to tho wall, nn offootlvo coll npHO angle enn bo 
calculated (nrotnn h/b>. Thin wan done nnd tho rosultlng plot of collnpRo angle 
vr passage blockngo Ir flhown In figure 52. Abbott* r roflidtR arc nhown hr dntn 
polntR nnd tho original cavity model collnpRo function In ahovVn nn n dashed lino. 
Doth tho Abbott dntn nnd our collnpRo function ngroo In that low onvlty hclghtR (or 
RtopR) hnvo high collnpRo nngloR (very nbrupt collnpRo) nnd high envltlos hnvo low 
collnpRo nngloR (vory gradual collnpRo). There Is n rnpld transition from low to 
high collapse angle ns tho cavity height Incronsos, Tho qunntltntlvo lovols of the 
Abbott data and our own collapse function arc dlfforent and thoro arc sovernl 
possible explanations for this. The step of Abbott Is not exactly similar to the 
cavity In that tho stop represents a solid flow boundary. Abbott's data are also 
representative of fully developed flow to which the step Is perpendicular nnd which 
Is confined by sidewalls. In tho Inducer the cavities may collapse on the uncovered 
part of the blade where the flow Is not confined and the flow will contain velocity 
and pressure gradients. These flow differences would be expected to result In a 
quantitative difference In the collapse angle functions. 
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Figure 49. Comparison of Cavity Length-to-Maxlmum DF 91080 
Height Predictions With Exact Predictions; 
ideal Hydrofoil Cascade, l = 0.78 rad (16 deg) 
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Figure 60. Comparison of Cavity Length to Maximum DF 91081 
Height Predictions With Exact Predictions, 

Ideal Inducer Cascade 
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Figure 61. Flow Separation Downstream of a Back- DF 91079 
ward Facing Stop (Data From Abbott and 
Kline) 
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Figure 52. Flow Separation Data Compared With DF 91152 

Empirical Cavity Coliapse function 

The general agreement of the cavity collapse angle function with the Abbott 
data was considered to substantiate the form of the cavity model collapse angle 
function. Accordingly, the original angle function was adjusted to bring overall 
cavity length predictions Into agreement with the test data. The final function Is 
shown as a solid line In figure 52. Note that the use of overall cavity length 
measured data to de+’ine the function effectively compensates for the Inaccuracy In 
length to maximum cavity height predictions that were discussed In the previous 
subsection. The error in profile area that Is accepted through this approach can 
be appreciated by considering figure 53, which shows cavity profiles In an Ideal 
Inducer cascade. The solid profile Is the one that would be predicted by the cavity 
model and the dashed line Is the Stripling profile tc maximum height along with 
the linear collapse to the same point as that of the cavity model. The area error 
Is approximately 13% for the small and 12% for the large cavity, values which are 
considered quite satisfactory. 

5. 1. 2. 4 Cavity Overall Length 

Predicted overall cavity lengths, generated using the collapse model of 
figure 52, are compared with Inducer test data from (1) In figures 54 and 55. 
Figure 54 Is for the Inducer tip measuremeni streamline and figure 55 Is for the 
midspan streamline. The available test data are limited and the cavity model 
cannot be fully substantiated, particularly with regard to the length variation 
with flow coefficient. A degree of confidence In the model's prediction as a 
function of flow coefficient is gained, however, when the correlation for both the 
tip and midspan streamlines Is considered. There Is a difference In Incidence 
angle at these two streamlines because of the difference In blado angle and wheel 
speed. The Incidence and blado angle differences would be predicted to have a 
significant effect on the cavity profiles, as would a difference In flow coefficient. 
A tabulation of blade and flow angles at those two streamlines Is given In table 8. 
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Figure 54. Comparison of Cavity Length Predictions DF 91077 
With Inducer Test Data, Tip Streamline 


* 

know* 

0 txow 
00,0*0 
— . Predicted 


o.oi aOJ 


iu txa tu 

CAVITATION NUMBKH. k 


Figure 65, Comparison of Cavity Length Predictions DF 91078 
With inducer Tost Data; Midspan Streamline 
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Figure 50. Comparison of Cavity Length Predictions 1)F 89779 
With Cascade Test Data 
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Figure 57. Predicted Suction Surface Cavity Volume, DF 90999 
Radial Leading Edge Inducer 

Sweepback tends to reduce cavity length as shown In the measured data plot 
of figure 58. It can be assumed that volume Is proportional to length so that 
volume also decreases with sweepback. This relationship of cavity volume to 
sweepback was used to evaluate Instability data for the swept Inducer Instead of 
a more rigorous regeneration of swept Inducer cavity volumes and stability 
predictions because of convergence problems which were being experienced with 
the Inducer Internal flow computer program. The computer program was being 
refined to Include swept Inducers concurrently In a related program (26). This 
approach was adequate for our purposes where an Instability prediction technique 
was to be validated and the cavity volume predictions for swept Indtcers will be 
obtainable when the work of (26) Is completed. 

5.2 CASCADE TUNNEL INSTABILITY 

5.2.1 Analytic Model 

A dimensional schematic of the cascade tunnel Is shown In figure 59. A 
linear dynamic model of the facility was prepared to Investigate the character- 
istics of cascade flow that could cause an Instability. The model was based on 
the assumptions that the flow past each cascade blade Is Identical, causing the 
cavities to oscillate In unison; and blade-to-blade pressure and flow distortion 
and cavitation propagation along the cascade do not occur. No propagation was 
evident during the tests. 
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A simplified model of the complete facility wa3 prepared so that the essen- 
tlal features of an Instability could be Investigated, Referring to the schematic 
of figure 59, It is apparent that the plenum tank, being large and containing a 
free surface, could be considered a point of constant pressure (Infinite compliance). 
If some plane In the nozzle could also be considered a point of constant pressure, 
these two points become boundaries that Isolate the nozzle and cascade from the 
remainder of the facility. Only that portion of the facility from the chosen nozzle 
plane through the cascade to the plenum tank requires modeling. The validity of 
treating a plane In the nozzle as a constant pressure boundary was later sub- 
stantiated by the test data. 



Figure 58. Measured Effect of Sweep Angle on Inducer DF 90998 
Cavity Length 




A schematic of that portion of tho cascade facility chosen for modeling Is 
shown In figure 60. The model consists of? 

1. A constant pressure boundary at static pressure (p n ) across 
which a fluid of constant density (p) Is flowing. 

2. A rigid nozzle length (/[) and Inlet and oxlt cross sectional 
areas (A n and Aj) measurod normal to the flow direction. 

Flow through the nozzle (W[) is considered Incompressible, 
Motionless, and one-dlmcnslonal at each cross soctlon. 

3. A cascade and tailboard section with a total— to-total pressure 
loss (D) that Is dependent on the Instantaneous value of Inlet 
static pressure (pi), inlet flowrate (Wj), and discharge flow- 
rate (W d ). Pressure drop Is measured from the cascade 
leading edge plane to the tank pressure (pj.). 

4. A cavitation cavity on each cascade hydrofoil suction surface 

that results In a total cavity volume (V c ) that Is dependent on 
the Instantaneous values of Inlet static pressure (p { ) and '"let 
flowrate (Wi). 1 

5. A large plenum tank at constant pressure (pj). 

The model neglects the boundary layer bleed system, which has the effect of an 
upstream compliance, and neglects fluid compressibility and wall expansion In 
• the nozzle, which would require that the nozzle be broken into several calculation 
sections (or ’lumps”) for accurate modeling. However, It Is unlikely that either 
of the neglected effects would make a difference In determining whether flow 
through the cascade Is stable or unstable. 



Figure 60. Cascade Instability Model FD 61195 
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For small amplitude oscillations, dynamic bohavlor of the model can be 
described, and the onset of an Instability can be predicted, by troatlng linear 
perturbations about the steady-state operating point. A sot of linear algobrnles 
and differential equations results. Lotting (A) denote a linear perturbation and 
(•) denote the derivative with respect to time (d/dt), the equations (general and 
linearized), which represent the model, ares 

Nozzle: 



Cavity Volume: 


V c = function (p t , W t ) 


dV 


dV 


AV ° S5 ap l Ap l + dW° AW 


l 


Continuity: 

w, = w d 



AW, = AW d - „AV 0 
Cascade Pressure Drop: 


p l + 


W 


2pAf 


W. 


PA 


-=P t +D 

ad 


(2a) 


(3a) 


(4a) 


D = function (p t , W t , W d ) 

Ap l + Jw[ AW l + ^ AW d ( 6 a) 

The coefficients In each of the above linearized equations are taken to be constants 
evaluated at the steady-state operating point, where Wt » Wd = W, the steady-state 
flowrate. Substituting equation (2a) Into (3a) and (4a) Into (5a), and taking the 
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Laplaco transform (S » Laplaco variable) of each equation, a sot of equations Is 
obtained In which some convenient definitions can be made 



The following definitions are made: 


Inlet Resistance: 
W 


r i-p 


a 2 " .2 
A. A 
L l n J 


Inlet Inertance: 


_ r~i dx 

i J x 


Pressure Compliance: 

dW c 

C = -p — 2. 

P K dP t 

Flow Compliance: 
dV c 

c w = ^ 


aw, 


i 


Pressure Gr in: 


G 

P dP t 


Flow Gain: 


r W 
w ' pA 2 


an 

aw, 


Cascade Internal Resistance: 


R 


dD 
c w; 


( 6 ) 
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Substituting thoso definitions Into equations (lb, 3b, and 6b) allows the character- 
istic equation to bo obtained In the forms 


IV'lV 8 + (R c<¥p - C w> + G p L l ,S + f R c + «1 G P - °wl - 0 


(7) 


An operating point will be stable If all roots of the characteristic oquatlon havo 
negative real parts, Indicating a tendency of the flow to return to an equilibrium 
state If momentarily disturbed. Otherwise the operating point Is unstable. Be- 
cause this characteristic equation Is of the second order the roots will have 
negative real parts (Indicating a stable operating point) only If all three coef- 
ficients In brackets have the same sign. Conversely, the operating point will 
be unstable If any coefficient vanishes or takes on a sign different from the other 
two. Therefore, to determine stability, It was necessary to evaluate each term 
making up the coefficients of equation (7). These terms were evaluated at the 
steady-state operating point from the definitions of equation (6). 


Inlet Resistance and Inertance (Rj and L t ): 
If Inlet resistance Is written In terms of the 


q 


1 


W 2 

2pAf 


Inlet velocity head (q^) where: 


The expression for resistance Is : 


R 


-fa 

l “ w 



Since both terms are positive It Is apparent that Rj Is always positive. 
When inlet Inertance Is defined In terms of a ’’mean nozzle area” (J) such 


L 


l 



h 

A 


It Is also apparent that Lq Is a positive number at all operating points. 
Pressure and Plow Compliance (C p and C w ): 

According to the cascade flow model of Stripling and Acosta (19) for fixed 
geometry and fixed Incidence, the maximum cavity height (h) and length (c) 
from the blade leading edge to the point of maximum cavity height (figure 46) 
are functions only of cavitation number (k) where: 

, Pi - p v 

k qj — (P v = vapor pressure) (8) 


In (19), both h and c are shown to Increase as k Is lowered and to decrease 
as k Is raised. The experimental measurements of cavity height and total 
length show the same trend. It Is reasonable to expect that total cavity 
volume (V c ) will also be a function only of cavitation number, for fixed 
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geometry and Incidence, and will follow a trond Identical to the longth and 
height as shown In figure 61. If wo assume that tho cascado Incidence Is 
fixed at any operating point, whothcr stoady or unsteady, then tho torms 
Cp and C w can be calculated (using equation (8)) ns; 


dv 

P K dp 




l ‘‘l 


dk 


c • -p— = 2£!s 
w P W dk 


Since dV c /dk Is always negative, C p Is always positive and C w Is always 
negative, * 




CAVITATION NUMBER, k 
Figure 61. Cascade Cavity Relationships 



FD 62908 


Pressure and Flow Gain (G p and G w ): 

The cascade total pressure drop (D) at any cavttatlng steady-state operating 
point can be written as the product of the noncavltatlng pressure drop (D nc ), 
at the same flowrate but high Inlet pressure, and a dimensionless pressure 
drop ratio (x ). 


This equation serves as the definition of X and as the means for calculating 
Its value at any operating point. The noncavltatlng pressure drop (D nc ) 
will be a function only of flowrate (W) for fixed geometry and Incidence. 

We assume that X is a function only of cavitation number (k). Then at any 
steady-state operating point: 

D(k,w)= X(k)D nc (W) 

To adapt this equation to the cascade Instability model, we assume that the 
same relationship holds under nonsteady conditions, with X always equal to 
the steady-state value of x that would exist at the Instantaneous cavitation 
number, and D nc always equal to the steady-state value of D nc that would 


exist at the Instantaneous discharge flowrate. Then at any unstoady opera- 
ting points 


n(k, W d ) « A(k) n nc (W cl ) 


dD „ n a A _ J) nc dA 
dPj ' nc c>Pj q t dk 

|a oD ax 

3Wj nc dW t W dk 

Defining a noncavltatlng prossuro loss coefficient (C nc )s 


D nc ® C nc q l 
G P = l -a^ =1 - C nodk 


clX 


G « 
w 


_SL . .12. = r, + k c ax.] 

pA 2 aw, W |^ 1K '-nc dk 


G« and G w can be positive or negative depending upon the sign and magnitude 
of (dA/dk), 

Cascade Internal Resistance (R c ): 
dD 

■ = A 

'd 


R, 


dD . 

aw. = A_ aw 


nc _ x ^ (C nc q lL 2xC 

X aw 2 nc W 


i 


R c Is positive at all operating points. 

The terms making up the coefficients of the characteristic equation (7) can be 
summarized as follows: 


R 


2q 


l 


l W 


l 

1 

. n J 


h -i*. 

1 A 

c .rfifTa 

P q t die 
c 

w W dk 


>0 

>0 

>0 


<0 


( 9 ) 


G„ 


G 


W 


H. 


" 1 " c no TiF 

O 

AV 

3r ItU c 

W [ l ' k L ne dfcj 

AV 

o 

q l 

nc W 

>0 


depend upon tho sign of 

dX 

dk 


Tho first coefficient In tho characteristic equation (7), H 0 L(Cp, Is therefore 
positive at all operating points and the other two coefficients, [H C (R[(.' - C’ w ) + g d Li) 
and [He + UiG p - G w ] must also bo positive for an operating point to bo stablo. Thero- 
fore, an operating point will bo unstablo If: 


K 0 (Hi c p - C w ) + Gpl^ < 0 


(10a) 


or If: 


U c +R t G p - G w < 0 (Ha) 

Referring to equations (9), (10a) and (11a), It Is apparent that the only way an 
Instability can occur Is for G p to be negative In equation (10a) or for R»G D - G w 
to be negative In equation (11a), since L[, R c and RiC d - C w are always positive. 
Therefore, the model Indicates that cascade Instabilities result from the gains 
(Gp and G w ), which are related to the term dx/dk, the slope of the pressure 
drop ratio curve. Stable and unstable forms of this curve are shown In figure 62, 




Figure 62. Cascade Pressure Drop Ratio Curves FD 62909 

(Cavltatlng/Noncavltatlng) 

In the stable characteristic, the curve Is flat at high cavitation numbers 
and turns upward smoothly as cavitation number Is lowered, Indicating that the 
cascade becomes a greater restriction to flow as cavitation Increases. In the 
unstable characteristic, tho curve Is alEo flat at high cavitation numbers, but 
has an Inflection as cavttatton number Is lowered. The requirement that the 
loss curve have an Inflection for a cascade Instability to occur was also found 
by Acosta (6). The slope of a loss curve Inflection that would be required to 
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cause nn Instability won defined by substituting from equation (9) Into equations 
(10a) and U In). An operating point would bo unstable (ft 

dx > 2XA 

dk " c p A 2 

no * j A| 

or If? 



1 ho right Hi do of oquntion (10b) Is posltlvo nt nil operating points. In oquntion 
(lib) the ordor of magnitude of the terms Is: 

A r 'l 



c 


nc 


rV> 


i 


A ^ A n ~ 0 * 01 

so that the right side of equation (lib) Is also positive at all operating points. 

Therefore, the model predicts that an operating point can be unstable only 
If dX/dk Is positive and has a magnitude greater than that required by either 
equation (10b) or (lib). Referring to the pressure drop ratio curves of figure 62, 
the model predicts that a cascade with the stable characteristic will be stable at 
811 fP^tlng Points, while a cascade with the unstable characteristic will be 
unstable on that part of the Inflection having a positive slope, If the slope Is 
sufficiently steep. The magnitude of the slope required for an Instability can bo 
estimated by approximating the values of the terms. These wore estimated as 
follows : 

X ~ 1 


A f /A n ~ 0 * 01 
k ~ 0,4 

A ~ 600 cm 2 (93 In?) 

~ 1,6 m (5,3 ft) 

A t ~ 500 cm 2 (78 In?) 

dV c /dk ~ -2000 cm 3 (-122 In?) 
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The value of d\' ( ,/dk wan estimated by fitting a circular arc from (lie blade 
loading ed^o through (he maximum height point, defined by e and h In flpyurn 4f», 
and back to the blade to define a completely circular cavity. Then, id an 
Incidence of 0,2B rad (10 den), the total cavil, v volume wan calculated for 
k • • 0,4 and k : 0,404, and d\' ( ./dk wan approximated iri AYT/'Ak, Snbatltutlni', 
the above valuea Into eouidlon (10b) and (1 lb) we obtain the follow Inf, call mates 
for the alope re(|ulred to cauae an (notability. 


The required slopes arc approximately equal, ho that the Inatablllty model 
predletlonH can bo summurl/.od as staling that an operating point In the neighbor- 
hood of k 0,4 will be unstable If dX/dk has a positive slope greater than about 
1,0, We would therefore expect to find an Inflection In the experimental pres- 
sure drop ratio curves, similar to figure 02, at each of the four test points 
shown as unstable In figure 15. 

5.2,2 Comparison of Cascade Instability Predictions to Tost Results 

The measured cascade pressure drop ratios from Inlet total to the plenum 
free surface are plotted against cavitation number In figure 03, Estimated maxi- 
mum uncertainties In cavitation number (J 2' V ) and In pressure drop (d-lO'V) are 
sufficiently small to allow Identification of Inflections In the curves. In figure 03 
there Is an Inflection In the curve for 0.28 rad (10 deg) Incidence at a cavitation 
number of about 0.4, and In the curve for 0.35 rad (20 dog) Incidence at a cavita- 
tion number of about 0.35. The curve for 0.31 rad (18 deg) Incidence does not 
show an Inflection. This Is probably due to an Insufficient number of data points 
between test points No. 20 and 1.7, where an inflection would be expected. The 
positive slopes of the two Identified Inflections arc much greater than 1.0, so, 
according to the instability model prediction, we would expect to see instabilities 
at all operating points that lie on the positive slope portions of the Inflections, 
Therefore, at 0,28 rad (10 deg) Incidence, unstable operating points are pre- 
dicted to be No, 2, 3, 4, 6, 6, 9, 36, and 30, and at 0.35 rad (20 deg) incidence, 
No, 12 and 13. All other operating points at these two Incidences are predicted 
to be stable. Definite predictions cannot bo made for 0.31 rad (18 deg) Incidence 
since the Inflection cannot bo Identified; however, It probably occurs between test 
points No. 26 and 27, and wo would expect test points to the loft of No. 26 and to 
the right of No, 27 to be stable. 

Comparing the predicted unstable points with the visual observations of 
table 4, most of the predicted unstable points are listed as unsteady, very unsteady, 
or pulsating. Howovcr, No. 9 and 30 are listed as steady. Also, many of the 
predicted stable points are llBtcd os very unsteady. There Is, therefore, no 
apparent correlation betwoon predictions and visual observations. This might 
have been expected slnco there was also no correlation between observations 
and measured pressure amplitudes. 
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Figure 63, Cascade Pressure Drop Hntto, Cavltatlng DF 91001 
to Noncavltatlng 

The test data points (28, 36, 37, and 40) with the highest cascade Inlet 
dynamic pressure amplitudes (0.28 to 0.46 N/CM 2 , 0.4 to 0.7 pst) are shaded 
oti figure 63, Points 36 and 37 could fall on the positive slope portion of the 
0, 28 rad (16 deg) Incidence curve and might be expected to be unstable. Points 28 
and 40 are removed from a positive slope area, however, and the relatively high 
dynamic pressure amplitudes for these points cannot be explained by our Instability 
model. There are also many points on the positive slope portions of the curves ‘ 
that would be predicted to be unstable but which had relatively low dynamic 
pressure amplitudes. 


Since there was apparently little correlation between predicted unstable 
operating regions and the test data, an estimate was made of the dynamic pressure 
amplitude that might be expected at an unstable operating point. An "order or 
magnitude" estimate was sufficient for this purpose, and It was made using the 
linearized equations of the tunnel model. (It should be noted that a linearized 
prediction of pressure amplitude Is very approximate?, and that an analog solution 
would be required In a thorough anulysls. ) Since very small amplitudes were 
expected, the linearized equations for nozzle flow and cavitation number were 
used to form a relationship between changes In static pressure (Apt) and changes 
J.n cavitation number (Ak) ns follows? 


Ap, - -(n, + L,S) AW, 

4k = ^ *Pi -w Jw t 


Ap, 

aT 


©JR, + L,B) 


1 + 


w 


2q l 


(H, + L,S) 
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Ufling tho definitions of oquntlon (f» for (IR) and (L.) and with 
W ~ Ai(2 p cq) 


la 

Ak 





Assuming that cavitation number varies sinusoidally at the frequency of test 
point No. 36 (20 Hz), then using the same approximation previously used for 
an Incidence of 0.28 rad (16 deg), and additionally (qi *0.2 N/cm 2 ), and 
(A[ « 0. 05 m^), and (p« 1000 kg/m 3 ), an estimate of 0. 2 N/cm 2 was obtained 
for the magnitude of (Apj/Ak). 

At an unstable operating point (No. 36 for example) we could expect 
cavitation number to vary In figure 63 so that the instantaneous operating point 
oscillates from one negative slope portion of the curve, where operation Is stable, 
through the positive slope portion, where operation Is unstable, and over to the 
other negative slope portion, where operation Is again stable. This would require 
cavitation number to vary from 0. 39 to 0. 41. For this variation In cavitation 
number, Ap( - 2(0.2) (0.41 - 0.39) = 0.008 N/cm 2 (peak-to-peak) (0.012 psl). This 
value Is approximately two orders of magnitude lower than any recorded dynamic 
pressure data. 

An estimate of expected cavity length variation was made by assuming that 
total cavity length was twice the length to maximum height (c) predicted by 
Stripling and Acosta (19). For swings In cavitation number between 0. 39 and 
0.41, total cavity length would vary from 44% chord to 42% chord, a movement 
that would be difficult to detect and that was considerably exceeded for all 
cavltatlng test points. 

Our conclusion concerning the cascade tunnel Instability correlation is that 
any Instabilities that may have occurred were obscured by normal cavitation 
oscillations. Cascade tunnel test results are therefore Inconclusive In regard to 
substantiation of the model as a means for predicting Instabilities. An analog model 
of the complete system would allow more accurate predictions of unstable frequency 
and amplitude, but the complexity of the system and the erratic, nonperiodic nature 
of the test data make It doubtful that such a model could be correlated with the data. 
Further modeling of the cascade tunnel was, therefore, not considered justified. 
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5. 3 INDUCER INSTABILITIES 
5. 3. 1 Analytic Model 

A drawing of the Inducer test loop, showing significant features and dimen- 
sions, Is shown In figure 64 and a schematic of the dynamic model that was pre- 
pared to represent the system Is shown In figure 66. The circled whole numbers 
refer to static and dynamic pressure measurement stations. Decimal numbers 
refer to calculation stations between measurement stations. Pressure drops 
around the loop are represented by resistances (R), fluid Inertia by Inertances 
(L), and fluid and pipe wall elasticity by compliances (C), using the standard 
lumped parameter modeling technique for hydraulic lines. Air pockets were 
treated as Increases In local compliance and, since a linear analysis was used, 
the resistances were linearized about the steady operating point. 

Frequency response analyses of the dead-end lines showed that the proper 
relationship of pressure to flowrate at the junctions could be obtained at frequencies 
up to about 30 Hz by treating the dead-end lines as simple compliances. The first 
closed-end resonance frequency of the longest section of the line In the Inner loop 
treated as a lumped parameter is 135 Hz, as determined from the one-dimensional 
wave equations. Generally, experience with dynamic models has shown that the 
lumped parameter modeling technique is accurate at f ?equencies up to 10% to 20% 
of the first resonance. The model would therefore accurately follow a sine wave 
input up to about 20 Hz, and become progress! vly less accurate at higher frequencies. 
This accuracy was considered sufficient to model the instability, since the funda- 
mental oscillation was found to be less than 20 Hz. 

Equations were written to define the "through flow" or branch portions of the 
test loop in terms of the measured change in static pressure between two points, 
as listed in equations (12) below. Compliance at the node locations indicated in 
figure 65 was accounted for by relating it to the difference between entering and 
leaving flow, as listed in equations (13), and the Inducer interfaces with the test 
loop system were defined through the inlet and discharge mass flowrates (Wg) 
and (Wg) and the inlet and discharge static pressures (pg and pio), as shown in 
equations (14) and (15). All equations represent linear perturbations of pressure 
and flow about the steady-state operating point, the conventional A having been 
omitted for clarity. S = Laplace variable representing d/dt. Other symbols 
are defined in the Appendix. 


















p io - P 1 - »10 W 10 - L 10 ^10 ' 0 

Pl-P 2 -W L l SW l“° 

p 2 ” P 2. 5 ” R 1 ' L 2 SW 2 " 0 
P 2 .5-P4- R 2.8 W 2. 6 - L 2.6 SW 2.5“ 0 

P 4 - P 5 - H 4 W, - L, SW 4 = 0 (12) 

P 5 " P 5. 5 " R 5 W 5 ” L 5 SW G 0 

p 5. 5 " p 7 " ‘li.S^t.S 511 

P 7 - p 8 -B 7 W 7 -L 7 SW 7 =0 

p 8- p 9- R 8 W 8- L 8 SW 8 =0 
W 9 - W 10 - C 10 ^lO ■ 0 

W 10 ” W 1 " * C 1 + C l. l> “ 0 
W 1 - W 2 - C 2 ^2 =l) 

W 2' W 2.6" C 2.5 Sp 2.5'° 


W 2. 5 ‘ W 4 “ C 4 Sp 4 = 0 

(13) 

W 4- W 6-< C 6 +C 5.X )Sp 6= 0 


W 5- W 6.5- C 6.5 Sp 5 .8=° 


W 6.6 ” W 7 " C 7 ^ =0 

. 


\ 

W 8 - W 9 - C w SW 8 ‘ % + C 9> Sp 9 = 0 

(14) 

p 10- G p p 9-°w W 8 + R p W 9“» 

(16) 


These equations constitute the Inducer system analytic model. 
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The values that wore calculated for (I,) and (C) from the known test 'oop 
geometry, tho calculated pipe compliance, and the bulk modulus of water at the 
moasu rod temperatures are tabulated In table 9, Resistances (R) at tho tost 
flow coefficients, which wore calculated from tho moasurod static pressures, 
are also tabulated. Discharge valve resistance was a variable (since Inducer 
head rise varied with cavitation number), and Its resistance Is plottod against 
cavitation number In flguro (iO for tho radial leading edge Inducor. 

The Inducer portion of the system dynamic model Is Illustrated In tablo 10. 
The variables that doflned the Inducor In equations (14) and (15) were obtained, 
as Indicated In tho table, from the moasurod head vs flow map (flguro 24), tho 
measured head vs cavitation number map (figure 27), and tho predicted cavity 
volume map (figure 57). Tho cavity volurno map of flguro 57 represents blade 
surface cavitation only, and, to account for tip vortex cavitation, It would be 
necessary to plot the sum of both forms of cavitation. The derivation of the 
equations that represent the Inducer are discussed In the following paragraphs. 


Table 9. Values of Inducer Loop Inertance, 
Compliance, and Resistance 


Position 

(Subscript) 

Inertance, 

L 

N/em 1 2 

Compliance, 

C 

kg/s 

Resistance, 

R 

N/cm“ 

kg/s 

kg/s 

(N/cm 2 )/s 

4> = 0. 070 

<t> = 0. 084 

<t> - 0.090 

1 

8.36 x 10~ 3 

0. 54 x 10" 4 

(1) 

(1) 

(1) 

1.1 

(2) 

5. 56 x 10" 4 

(2) 

(2) 

(2) 

2 

11.20 x 10~ 3 

0. 76 x 10" 4 

4.66 x 10“ 2 

5. 64 x 10‘ 2 

6.10 x 10" 2 

2.5 

15. 00 x 10“ 3 

1. 01 x 10~ 4 

11.58 x 10" 2 

13.91 x 10~ 2 

15.01 x 10” 2 

4 

16. 00 x 10" 3 

1.35 x 10~ 4 

-3. 01 x 10“ 2 

-3. 51 x 10 -2 

-3. 76 x 10~ 2 

5 

6.40 x 10" 3 

1.45 x 10" 4 

2.71 x 10~ 2 

3.01 x 10‘ 2 

3.52 x 10“ 2 

5.1 

(2) 

19.30 x 10“ 4 

(2) 

(2) 

(2) 

5.5 

4. 10 x 10’ 3 

3.80 x 10' 4 

(2) 

(2) 

(2) 

7 

5. 49 x 10" 3 

19. 90 x 10“ 4 

1.07 x 10" 2 

1.29 x !0“ 2 

1.41 x 10“ 2 

8 

2. 71 x 10 -3 

3. 04 x 10~ 4 

0. 54 x 10~ 2 

0. 65 x 10" 2 

0.71 x 10~ 2 

9 

(2) , 

3.32 x 10" 4 

(2) 

(2) 

(2) 

10 

3. 87 x 10“ 3 

3.21 x 10~ 4 

5. 11 x 10“ 2 

6. 01 x 10‘ 2 

6.56 x 10 -2 


(1) Variable. (Seo figure G6. ) 

(2) Not Used In Model. 
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Figure 66. Discharge Throttling Valve L.-near DF 91007 

Resistance, Radial Inducer 


Under steady-state conditions. Inducer performance can be written 
(assuming constant speed): 


X = f^Pg.W) (16) 

P rnc = *2 <"> 

Under dynamic conditions. Inlet flowrate (W 8 ) can differ from discharge flowrate 
(W9) because cavity volume can change. Inducer pressure rise becomes a function 
of both flowrates. It Is assumed that the above relationship for Inducer pressure 
rise Is valid under dynamic conditions If the steady-state flowrate (W) Is replaced 
by Inlet flowrate (W 8 ) In "fi," and by discharge flowrate (W 9 ) In "f2. " Then 
under dynamic conditions Inducer performance can be written: 


x - flaw's) (17) 

P rnc = f 2<W 9 ) 
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Table 10, Inducer Representation In Dynamic Model 



Performance 
Continuity 
Cavity Volume 

Pressure Gain 
Flow Gain 
Pump Resistance 
Pressure Compliance 
Flow Compliance 


^10 = G p^9 + G w AW 8‘ U p AW 9 
AW g - AW (J = - PAV c 


-*av c = c p ap 9 + c w aw 8 


d p 10 

G = = 1 + 

p a p 9 


2± 


1 +4> 


bit 




Head 

Falloff 

Map 


Head 

Flow 

Map 


Cavity 

Volume 

Map 


In these equations, "fj and f 2 " still denote the her.d falloff and head vs flow maps 
obtained from steady-state tests or predictions. For linear perturbations about the 
steady-state operating point. Inducer discharge static pressure (Pio) can be written: 


AP 10 - A p 9 + AP r ■ -m 8 


dP r 

A p 

r OPq 


APo 


+ Aq fl ■ Aq 9 

aw Q Aw s + aw. 


AW n 


(18) 

( 19 ) 
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Tho symbol (A) denotes a llnoar perturbation about tho steady-state operating 
point. Inlet and discharge velocity heads are denoted by qg and q 9 , respectively, 
Kach coofflclont In equation (19) Is considered to bo a constant, evaluated at tho 
steady-state operating point. From tho relationships of equation (17), tho coef- 
ficients In equation (19) can bo evaluated as? 



P 

_&L 

df, 

r- f .. . -i 


rnc 

dPc, 

2 dp 

3l> r 


dx 

f 3f l 

aw 8 • 

= p 

rnc 

w 8 

sf 2W 

&>„ 

dP 



5% ' 

v rnc 

^7 “ 

f llW 


(20) 


The velocity heads can be written: 


«8" 


W 2 

8 


2 PA 


l 


<*9 = 




2 PA 


Aq g = 


d 

W 


PA 


1 


AW C 


( 21 ) 


Aq 9 = 


W 


PA 


2 AW 9 


Substituting equation (20) Into equation (19) and equations (19) and (21) Into 
equation (18), Inducer performance becomes: 
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The throo coefficients In oquutlon (22) can bo Intel pro tod nn a "pressure 
gain" (Gp), n "flow gain" (G w ), and a "pump Internal resistance" (It,,), no that 
equation (22) can bo written! 1 


APlO " V P 0 ' G w AW fi " V W 0 


(23) 


The throo coefficients could bo ovaluatod dlroctly from tho bond fa 1 1 off and 
bond vs flow maps and known operating point If tho maps wore available In torms 
of pressure and flowrate, Since tho maps arc available In terms of cavitation 
number (k), flow coefficient (£), and noncnvltatlng hoad coefficient ($), It Is 
convonlont to rodoflne tho functions "fj and f 2 " In those torms, 


\ E f^(k,(j>) 

£ * f 2 (5) 

Use Is made of tho following standard relationships! 


(24) 


k = 


p 9~Py 

q' 


$ = 


U 


t 


v 51 

V l - PA, 


(25) 


q’ 

P 


— (U 2 + V 2 ) 

2' t V 

- 


1/2 


rnc 

The coefficients (G p , G w and R p ) In equation (23) can then be obtained In 
the equivalent form shown In table 10. Values for these coefficients that were 
calculated for the radial leading edge Inducer are given In figures 67, 68, and 
69* 

Continuity at tho Inducer Involves changes In the total volume (V c ) of 
cavitation present. For Unear perturbations, the liquid being stored In tho 
Inducer Is: 


* W 8- AW 9“dt<-* AV c) 


(26) 


Tho Instantanoous cavity volumo Is assumed to be a function (fo) of Inlet 
static pressure and Inlet flowrate, where "f 3 " denotes the steady-stato cavity 
volumo map. 


V c = f 3<P9’ W 8> 


(27) 
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For Unonr perturbations: 


av 


av. 


AV » - 1 


Tho liquid being stored In the Inducer hooomnst 

av , av 

- pAV c- p ept J: 


<8«) 


(UK) 


l ho two coefficients In equation (2!)) can bo Interpreted aw a pressure com 
pllanoo (( p) and a flow compliance (C w ), ho that o(|uatlon (2D) can bo written: 


" ^ V c C p 


^9 , C w^ W M 


(tin) 


Tho two coefficients, Cp and C^, could bo evaluated directly from the 
cavity volume map If tho map wore available In terms of pressure and flowrate. 
Since the map Is available In terms of cavitation number a.d flow coefficient, 

It Is convenient to rodoflno tho function "f a " In those terms: 

V c af g (R,5) (81) 

The coefficients Cp and C w in equation (30) cun then be obtained In tho 
equivalent form shown In tublo 10. Values that were calculated for tho radial 
loading edge Inducer are given In flguros 70 and 71. 


(0.02? J) 1 



Figure 70, Pressure Compliance, Radial Leading DF 91010 

Edge Inducer 
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Inducer 


Equations (23), (20), and (30) completely dol’lno tho Inducer dynamics and 
complete the model description. Predictions were generated by simultaneous 
solution of tho Inducer equations and equations (12) and (13) that describe tho 
tost loop hydraulic system. 

5.3.2 Model Predictions 

Model predictions were generated In tho form of a damping ratio ({*) and 
an undamped natural frequency (u» that corresponded to each root In the character- 
istic equation of the water loop. Graphical definitions of f and u> are given In 
flgut c 72. Tho roots Indicate tho manner In which the water loop would respond 
to a disturbance. If "y(t)" represents an arbitrary parameter In the water loop, 
such as Inducer discharge pressure, each pair of complex roots Indicates an 
oscillatory response of the form: 

y(t) - e sin (to/ 1 - t 2 t) 

and each real root indicates a ramping response of the form: 
y(t) ’• e “ rwt 

The overall response is the sum of the response from each root. Tho actual form 
of the disturbance (sine wave, pulse, etc. ) would add additional terms to the over- 
all response; however, the roots Indicate the response that Is characteristic of 
the water loop and determine whether the water loop would respond in a stable or 
unstable manner. Therefore, to analyze the stability of the water loop It Is only 
necessary to examine the roots of the characteristic equation for evidence of 
unstable response. The disturbance required to Initiate unstable response Is 
always presort In the form of normal vibration and need not bo considered. 
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Positive damping ratio (f) indicates a decaying oscillatory amplitude or a 
decaying ramp, which is stable behavior. When damping ratio is negative, the 
water loop is predicted to be unstable; and when it is positive, the loop is pre- 
dicted to be stable. The form of the unstable behavior is indicated by whether 
the damping ratio is obtained from a complex root pair (oscillatory behavior! or 
from a real root (ramping behavior). 

The convention used is that damping ratio always jies between -1. 0 < f < 1. 0. 
For a complex root pair, damping ratio can be interpreted as the ratio of actual 
damping to critical damping of a second order system, such as a spring-mass- 
damper sjstem, whore critical damping is the value at which oscillations would 
cease. For a real root, damping ratio is either r-1. 0, indicating whether the root 
lies to the left or right of the origin (figure 72). 

Actual osclllatory _frcq uency for a complex root is given by the damped 
natural frequency fcn/ 1-f 2 ); however, since the magnitude of damping ratio 
for the complex roots is generally very small, the frequency can be approximated 
by the undamped natural frequency ((*)), For a real root, natural frequency does 

not imply oscillatory behavior, but locates the distance the root lies from the 
origin (figure 72). 


Since the equations usod to describe the water loop are linear, amplitudes 
reached at unstable operating points cannot bo predicted. A linear analysts can 
only determine whether a given operating point will be stable or unstable, and, 

If unstable, Indicate whether the Initial bohavlor will be an oscillation or a ramp. 
There is no direct relationship between damping ratio obtained from a linear 
analysis and oscillation amplitude reached by the actual nonlinear system. 
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Howovor, a qualitative relationship usually exists In that a largo negative damping 
ratio Indicates that oscillation amplitude would grow much more rapidly than a 
small nogatlve damping ratio, Indicating a tendency toward a more severo 
Instability. It Is therefore reasonable to expect that damping ratio will be an 
indicator of relative oscillation amplitude. 

Actual oscillation amplitude can bo predicted by solution of the nonlinear 
equations describing water loop dynamic motion on an analog computer. Non- 
linearities arlso from pressure drops around the water loop and from the nonlinear 
relationship between cavity volume, pressure, and flowrate. A linear analysis was 
selected over the nonlinear approach because it allows greater flexibility in 
identifying destabilizing effects throug parametric studies and Is, therefore, 
a more powerful tool for evaluating mociu. validity. An additional nonlinear 
simulation would have complemented the analysis but was outside the scope of 
effort. 

Instability predictions were generated for the radial leading edge inducer 
and later qualitatively extended to the swept inducer during the data correlation 
effort. The only significant difference between the inducer systems was the 
smaller cavitation volume of the swept inducer. 

Noncavitatlng values of f and w are listed in table 11 in order of increasing 
frequency. Damping ratios are all positive and therefore correctly Indicate that 
the system is stable under noncavitatlng conditions. Roots 1 and 10 are real 
roots and represent ramping behavior of the water loop, while roots 2 through 9 
are complex roots and represent oscillatory behavior. Root 2 represents the 
basic low frequency, oscillatory behavior of the loop and was found to show the 
greatest change with cavitation number. Root 1 was found to indicate a runaway 
instability in the head breakdown region. 

Table 11. Noncavitatlng Values of Damping Ratio and Natural 
Frequency, Radial Leading Edge Inducer 


Root 

Number 

0 = 

0.070 

0 = 

0.084 

0 “ 

0.090 

Natural 

Frequency, 

Hz 

U) 

Damping 

Ratio 

r 

Natural 

Frequency, 

Hz 

U) 

Damping 

Ratio 

r 

Natural 

Frequency, 

Hz 

U) 

Damping 

Ratio 

r 

1 

4.5 

1. 0000 

3.6 


3.3 

1. 0000 

2 

44. 7 

0.0180 

44.7 

0.0152 

44.7 

0.0148 

3 

63.2 

0. 0525 

63.2 

0.0531 

63.2 

0. 0536 

4 

110 

0. 0209 

110 

0. 0125 

110 

0.0098 

5 

145 

0. 0309 

145 

0.0316 

145 

0.0321 

6 

172 

0.0008 

172 

0.0009 

172 

0.0010 


196 

0. 0140 

195 

0. 0094 

196 

0. 0079 


245 

0. 0102 

245 

0.0102 

245 

0. 0104 

9 

295 

0.0124 

295 

0. 0072 

295 

0.0055 

10 

1106 

1. 0000 

1100 


1097 

1.0000 
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Oscillatory type Instability predictions are shown in flguros 73, 74, and 75 
over a range of cavitation number from 0. 22 to 0. 04 for flow coefficients of 0, 070, 
0.084, and 0,090. Hoots 2, 3, G, and 8 hnvo nogatlvo damping ratios (an instability 
is predicted) over nearly the entlro range of cavitation number. Other oscillatory 
roots had positive damping ratios. Inducer hoadrlse was assumed to be unaffected 
by cavitation over this range, and the negative damping predictions are, thoroforo, 
representative of "continuity" mechanism instabilities. The measured shape of 
the head falloff curvo did not significantly affect predictions in this range of 
cavitation number. The air that was observed to form in the heat exchanger was 
accounted for in the analysis and is shown on the curves along with predictions 
made with no air. The noncavitating values of the roots are indicated on the 
right of the figures. 

Root 2 is the most Important because it represents the basic low frequency, 
oscillatory behavior of the system. Root 3 is only slightly higher in frequency, 
but its damping ratio does not show as substantial a change with cavitation number 
as that of root 2, Roots 6 and 8 are of such a high frequency that they are not 
considered significant to any instabilities. The root 2 damping ratio curve 
intercepts zero damping at a cavitation number above 0. 27. The curves show* a 
gradual drop in damping ratio with cavitation number to the point of air formation 
in the heat exchanger. At that point, there is an abrupt drop in damping ratio. 
Natural frequency also gradually drops with cavitation number to the point of air 
formation, at which point it drops abruptly. The natural frequencies and damping 
ratios generally decrease with flow coefficient, signifying that oscillatory 
frequencies would decrease and amplitudes become more severe as flow 
coefficient was reduced. 

Root 1, which represents a ramping type instability, also varied with 
cavitation number. Since this root is strongly dependent on the shape of the head 
falloff curve, the measured maps from figure 27 were used to generate predictions. 
Predictions are shown in figure 76 for the flow coefficients of 0. 070 and 0. 090; 

0. 084 was similar to 0.070 and was omitted for clarity. The predictions are 
shown in the form of damping ratio and natural frequency as a matter of convention. 
These terms are ambiguous for ramping-type instabilities, and they do not imply 
oscillatory behavior. A reduction in natural frequency toward zero Indicates a 
tendency toward transition from stable to unstable operation. The predictions 
therefore show that a ramping instability should occur for the 0. 090 flow 
coefficient at 1 = 0. 035. The other two flow coefficients are predicted to be 
stable, but they are tending to go unstable, as evidenced by the dropping natural 
frequency prediction. 

Model predictions are correlated with Inducer test data in the next para- 
graph, and the degree of correlation Ib evaluated through parametric treatment 
of the variables in paragraph 5. 3. 4. 

5. 3. 3 Correlation of Model Predictions With Test Results 

Model predictions for the lowest frequency oscillatory root (root 2) ar. 
compared with test data for the flow coefficients 0) of 0. 070, 0. 084, and 0. 090 
in figures 77, 78, and 79. Test data are from figure 39 and are the lowest 
measured frequency component. The measured range of cavitation number 
(0. 22 to 0. 03) and the measured head falloff curves were used in the model. 
Predictions are for the radial leading edge inducer, and test data are for both 
the swept and radial Inducers. The noncavitating predicted values of damping 
ratio (f) and natural frequency(tu) are Indicated on the right of the curves. 
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Figure 76. Linear Model Predictions of Inducer DF 91015 

Loop Ramping Characteristics, Radial 
Leading Edge Inducer 

The model predictions for the radial Inducer correlate favorably with the 
oscillatory test data in that: 

1. The negative damping ratio curve indicates that the loop is 
predicted to be unstable at all test points, and oscillations 
were measured at all test points. 

2. The predicted drop in damping ratio at k = 0. 06, when air 
formed in the heat exchanger, coincides with the measured 
sudden increase in oscillation amplitude. 

3. The predicted rise In damping ratio toward stable operation 
below 1c = 0. 04 coincides with the measured decrease in 
oscillation amplitude. 

4. The relative magnitude of damping ratio at the three flow 
coefficients when air formed in the heat exchanger agrees 
with the relative magnitude of measured amplitude. 

Predicted frequency is higher than measured and shows more change with 
flow coefficient than measured, but has the proper trend of decreasing frequency 
with decreasing cavitation number. 
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A runaway instability, ranking It Impossible to hold a stoody-stato operating 
point, occurred at a flow coefficient of 0,000 but not at the two lower flow coof- 
I'lclenlH. A real root of the loop characteristic equation (root 1), an Indicator of 
ramping behavior, predicted that a runaway Instability should occur at £ (), 0Q» 

but not at the two lower flow coefficients. The predictions were shown In figure 70. 
The figure shows an abrupt decrease In damping ratio at a cavitation number of 
a.oon for the n.000 flow coefficient, Indicating a transition to unstable ramping 
operation. The observed runaway Instability occurred at a cavitation number of 
approximately 0, o:u». 

The effect of leading edge sweepback on measured cavity length was shown 
In figure 68. The data arc limited, but they consistently show that the cavities 
become shorter as sweepback Is Increased, The effect of sweepback on system 
dynamics Is therefore to provide a smaller suction surface cavity at a given 
cavitation number. Swept Inducer system oscillatory characteristics would 
therefore be predicted to resemble radial Inducer characteristics If the two were 
compared with the swept Inducer operating at a cavitation number that was lower 
than the radial Inducer, such that cavity volume for the two Inducers was the same. 
This qualitative prediction Is consistent with the test data, which showed lowor 
amplitude pressure oscillations and higher frequencies at a given cavitation 
number for the swept Inducer. 

Prerotation of the Inducer Inlet flow can effoct stability through Its effect 
on the head flow map, the head folloff map, and the cavity volume map. Model 
parameters that were derived from these maps Include pump resistance (H p ), 
pressure gain (G p ), flow gain (G^), pressure compliance (C p ) and flow com- 
pliance (C w ). Since Hp, G p , ana G w were obtained from measured data, the 
effect of any prerotation that occurred in the tests Is already Included. C p and 
C w were obtained from the predicted cavity volume map, however, and do not 
Include prerotation effects directly. The cavity volume calculation does Include 
an empirical length adjustment, but the Inducer data used for the adjustment 
were Insufficient to empirically account for prerotation. Prerotation would tend 
to reduce Incidence angle and cause the cavitation cavities to be smaller than 
predicted, thereby reducing C p and C w . 

From figure 42, It can be seen that prerotation, as measured at a radius 
near the Inducer tip, was greater for the low flow coefficients and decreased 
as flow coefficient Increased. Actual cavity size would therefore tend to differ 
from the predicted "zero prerotation" cavities more at low flow coefficients, 
and there would be less spread In cavity volume with pressure and flow (C p and 
Cw) than predicted. The use of reduced C p and Cw values In the model equations 
would rosult In predictions that Bhow less spread In damping ratio and natural 
frequency with flow coefficient than the predictions of figures 77, 78, and 79. 

The new predictions would agree more closoly with test results that show very 
little spread In amplitude and frequency with flow coefficient at cavitation numbers 
(k) above 0. 06 (where prerotation was measured). Prerotation ceased below 
It = 0. 06, and a spread with flow coefficient In the test amplitudes and frequencies 
and In model predictions then occurred. However, below It = 0. 06 tho trend of 
test frequency with flow coefficient for the radial Inducer Is opposite to predictions, 
for reasons which are not presently understood. For the swept Inducer, the 
expected trend of Increasing frequency with Increasing flow coefficient occurred 
below It = 0. 06. Above It = 0. 06 the fundamental oscillation component for the 
swept Inducer was not sufficiently distinct to Identify Its frequency, so that a 
trend of frequency with flow coefficient could not be Identified. 
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Tlu> only significant nron In which llx* predictions do not correlate closely 
with (lie I csl did (i Ik the ni""til1ude of the oscillatory frequency, In nil effort to 
resolve the frequency differ* net*, the test procedure and facility were examined 
to Identify any unaccounted IV - system compliances that could have affected 
predicted frequencies, Itefoi being filled, the water loop was evacuated to 
1 N/env (It pslat, Previously deaerated water was then pumped Into the loop 
until loop pressure reached 17 N/enr C5 pslat. This pressure was maintained 
while bleed valves In high ptuMs of the loop were opened to expel trapped air, 

A review ol the bleed valve locations revealed that one area, in an unused loop 
ley;, could have had a trapped air pocket because the valve location was not at 
the (op of the pipe, The section of Pie loop In riuestlou Is shown In figure ho, 
Approximately ratlin env 1 (1..TI gnli of air til i N/em- (;i pslat could have been 
trapped by the venturi throat and compressed to <i:»n cm :i (n, :‘(l gab at 17 NAmi- 
(!i5 pslat dll rl lip the hi ceil process. These water levels are shown In the pipe 
cross section of figure si. The air pocket would remain trapped by the venturi 
tliront since loop pressure was not again lowered to 2. 1 N/cm- during tin* test. 
Hitch an nmounl ol air Is consistent with ail observed Increase In total loop fluid 
volume (as determined by a change In discharge accumulator surge tank level) 
of approximately 4785 env 1 (1. 0 gal) when loop pressure was lowered from 
110 to la N/eirm (OH to 15 psla). The trapped all* would account for 055 env* 

(0. 17 gal) of this volume change. 

The effect of the trapped all* on the system was calculated al a flow coef- 
ficient of n. 081 and a cavitation number of 0. 10. Trapped ah* volume under 
these conditions would be sno env* (0. M gal). Model predictions of damping 
ratio and natural frequency as functions of all* volume at the location shown In 
figure 80 are shown In figure 82. Air volume has a negligible effect on damping 
ratio, but It has a large effect on frequency. The calculated amount of alt* would 
lower the system natural frequency from 42 to 27 Use, ns compared with a 
measured frequency of 10 llz. As little ns 15% of the calculated amount of all* 
would lower the natural frequency from 42 to 54 Hz. It Is probable that nearly 
the total calculated amount of all* was trapped and that the all* accounts for part 
of the difference between predicted and measured frequencies. 



Figure 80. Location of Trapped Air Pocket FD G23G8 
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Figure 8 J, Pipe Section at Trapped Air Pocket FD 62369 



Figure 82* Effect of Trapped Air On Oscillatory DF 91005 

Predictions; Radiol Leading Edge 
Inducer; 0 = 0,084, li <= 0,16 
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Thus, there remains approximately 2900 em 3 (0,8 gal) of unue counted for dis- 
placement, Indicating the presence of additional compliances, which, If Included 
In the model, would certainly lower predicted natural frequencies, Tip clearance 
cavitation would form a part of this unaccounted for compliance; H i effect on 
damping ratio and frequency are shown In figure 8.‘l for the 0 0, 084, E - 0 , 1G 

condition. The effect of tip clearance cavitation would bo to Increase pressure 
compliance an amount that Is currently Indeterminate. An increase from the 
nominal value of 0.00007 (which represents predicted blade suction surface 
cavitation only) to 0, 084 would lower the natural frequency from 28 Hz to the 
measured 1G Hz. Intuitively, such an increase due to tip clearance cavitation 
alone seems high, but It should be remembered that there Is a substantial amount 
of unaccounted for compliance that would lower frequency regardless of Its source. 

+ 0 , 1 . 


Ortjtindl Ptcdklior 

Air t*oikcl In 

Itad-Und Vonturi tine (W0 tM4 gal) 


ItKttdalug Tip Vortex CavlUUon . 


E 20 

i - 

p Test irtqucno 

I 

0.0*101 0.000.’ 


0.0010 0.0030 0.0030 11 . 0100 " 

P.lMtmt' COM I’ll AMt, . Jl 

p nS 


Figure 83, Effect of Pressure Compliance on 

Oscillatory Predictions; Radial Loading 
Edge Inducer; 0 «= 0,084, h « 0,16 
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5. 3,4 Pnramotrle Effects of System Variables on Stability 

Several pnrnmotrlc studies woro e nductod during which the sonsltlvlty of 
Instability predictions to various systorn parameters was evaluated. The purposes 
of the studies woro to (1) provide quantitative Information concerning possible 
Inaccuraclos In tho predictions that may have rcsultod from Inaccuracies In 
parameter definition and (2) Identify, potentially "stabilizing" system changes. 

The results of thoso studies are discussed In tho following paragraphs, 

5. 3. 4. 1 Inlet and Discharge Line Incrtance and Resistance 

The effect of Inlet and discharge line Inertances (Ls and Ljo) on predicted 
damping ratio and natural frequency Is shown In figure 84 at a flow coefficient of 
0. 084 and a cavitation number of 0. 16. Increasing these inertances is equivalent 
to increasing the line lengths between stations 8-9 and stations 10-1 of figure 64. 
An inlet inertanco increase from nominal to five times nominal results in only 
slightly lower damping ratios and natural frequencies. The same increase in 
discharge Inertance has a negligible effect on the predictions. A similar lack 
of sensitivity was found for pump resistance (R p ), which Is proportional to slope 
of the head vs flow map, and for control valve resistance (Ri) over the same 
order of magnitude range, although these curves are not shown. These results 
are valid for our closed-loop system only. An open-loop system, such as a 
rocket engine feed system, was found to be sensitive to Inlet line length in 
another study (27). The above comments refer to oscillatory Instabilities. 

Ramping instabilities are sensitive to pump resistance (R p ) and control valve 
resistance (Rj) In that a decrease in these resistances reduces system resistance 
and head required to maintain flow. This reduction makes it easier for head 
available to jump above head required and initiate a runaway Instability. 
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8.3.4. 2 Cavitation Pressure Compliance (C p ) 

Pressure eomplnlnee Is a measure of the cavitation volume response to 
static pressure changes. The effect of pressure compliance on predictions was 
shown In figure 83 for $ = 0.084 and E - 0. 10. An Increase In compliance from 
the nominal 0, 0007 to 0,010 would cause system to be stable and would lower 
oscillatory frequency from 27 to 15 Hz. More than the indicated nominal pressure 
compliance was known to bo present In the inducer loop, the additional compliance 
being attributed to tip clearance cavitation, Futhor additions In pressure compliance 
could be obtained with an accumulator at the inducer inlet. A highly significant 
finding with respect to pressure compliance was that it alone cannot cause an 
Instability. If cavity volume changed only in response to changes in static pressure, 
the system would be stable. 

5. 3. 4. 3 Cavitation Flow Compliance (C w ) 

Flow compliance is a measure of the change in cavitation volume with 
flowrate. The effects of flow compliance at £ = 0. 084 and Tc = 0. 16 are shown 
in figure 85. A reduction of flow compliance from the nominal 1. 7 to 0. 25 ms 
would stabilize the system. Flow compliance has a negligible effect on frequency. 

Tip vortex cavitation would contribute some flow compliance above the nominal 
but the amount Is Indeterminate. Flow compliance is the predicted destabilizing 
parameter above head breakdown. If it is set to zero, the system la n 

to be stable regardless of the values of other parameters. 
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5. 3. 4, 4 Tnducer Pressure Gain (Gp) 

Pressure gain is a measui’e of the change In Inducer head rise with inlet 
prossure. The only area whore pressure gain was sufficient to significantly affect 
predictions was In the head falloff area. It was thorefore evaluated for flow 
coefficients of 0, 090 and 0 . 070 at a cavitation number of 0. 04. The sensitivity 
of oscillatory predictions to pressure gain In this area Is shown In figure 86 . The 
values of Gp for a 1c of 0, 04 and 0 , 031 arc indicated on the curves for reference. 
An increase In pressure gain Is equivalent to a reduction in cavitation number 
because the inducer head rise drops increasingly rapidly as cavitation number is 
reduced. As pressure gain is increased (cavitation number reduced), the model 
predicts a return toward stable operation. Frequency is predicted to decrease 
for 0 = 0 . 090 and to increase for 0 = 0 . 070. Pressure gain is the parameter that 
causes the predicted damping ratio in figures 77, 78, and 79 to rise toward stable 
operation below k = 0. 04, a prediction that compared favorably with measured 
data. 


Pressure gain has an overall stabilizing effect on ramping behavior of the 
test loop. 
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Figure 86. Effect of Pressure Gain on Oscillatory 

Predictions; Radial Leading Edge Tnducer: 
k = 0. 04 


DF 91018 


5. 3, 4. 5 Flow Gain (G w ) 

Flow gain Is a measure of the change In Inducer head rise with Inlet 
flowrate, The effects of flow gain on system oscillatory behavior are shown 
In figure 87 for flow coefficients of 0. 070 and 0. 090, and a cavitation number 
of 0. 04 and 0. 001 are shown for reference. The curves show that flow gain has 
a relatively small effect on oscillatory damping ratio and frequency. Flow gain 
has a strong Influence on ramping behavior, however, as shown In figure 88. 

The values of flow gain for the k = 0. 04 and 0. 001 points are Indicated. As 
cavitation number Is lowered from 0. 04 to 0. 001, flow gain decreases at 
0 - 0* 070 (and at 0. 084, jvhich is not shown) and increases at 0 ~ 0. 090. The 
Increase In flow gain at <t> = 0. 090 Is sufficient to cause the model to predict a 
transition from stable to unstable operation. The decrease In flow gain with 
cavitation number at the other flow coefficients is stabilizing. This predicted 
behavior correlates with radial inducer test data where a ramping-type instability 
was observed for <P = 0. 090, but not for other flow coefficients. Y 
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Figure 88 . Effect of Flow Gain on Ramping Pre- DF 91019 

dictions; Radial Leading Edge Inducer; 
k = 0.04 

In spite for the stabilizing effect of flow gain at 3> = 0. 070, a tendency toward 
a runaway instability below k = 0. 04 was predicted in figure 76. (The reduction in 
natural frequency Indicates that the real root damping ratio Is approaching a transi- 
tion from +1 . 0 to -1. (1 ) The reason for the predicted tendency toward instability 
was that inducer head rise decreased below k = 0. 04, which caused the control 
valve resistance (Ri) and the Inducer Internal resistance to decrease. The 
decrease in Rj and Rp is a destabilizing effect for the runaway instability since 
it reduces system resistance and head required to maintain flow. This reduction 
makes it easier for head available to jump above head required and initiate an 
instability. 

5. 3. 4 . 6 Inlet Line Compliance (C 7 ) 

The addition of inlet line compliance to the system model at a point some 
distance upstream of the Inducer (location 7) is comparable to changing the model 
from a closed to an open loop, wherein the inlet line becomes more similar to an 
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engine Inlet line. The effect of a high compliance (open loop) on damping ratio 
and natural frequency predictions Is shown In figure 89, 'Jhe solid line Is the 
original prediction (with no nlr In the heat exchanger), and the dashed lino Is for 
high Inlet line compliance. The prediction Indicates that the transition from 
stable to unstable for the two systems would occur at approximately the samo 
k (< 0,22), but that the open loop would probably oxpcrlonco groator amplitudes. 
The damping ratios and natural frequencies shown In figure 89 aro for the lowost 
frequency oscillatory root that showed significant change with cavitation number 
for each system. Unlike the closed-loop system, the next to lowest frequency 
root In the open-loop system (at high cavitation numbci’s) showed greatest change 
and became the lowest frequency root at low cavitation numbers. 
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Figure 89. Effect of Inlet Line Compliance on DF 91264 

Oscillatory Predictions; Radial Leading 
Edge Inducer, <p = 0.084 


5.4 SIMPLIFIED INDUCER SYSTEM MODEL 


A somewhat complex sot of equations (12, 13, — ) was used to describe the 
inducer test loop In paragraph 6.3.1. This complexity resulted In accurate re- 
sponse predictions but, unfortunately, It also obscured the significance of the major 
parameters affecting stability. A simplified system model Is described In this sec- 
tion to permit more apparent analytic interpretation of the roles played by pres- 
s urc gain (Gp), flow gain (G w ), pump resistance (R p ), pressure compliance (C n ), 
and flow compliance (C w ), The Influence of the hydraulic system parameters, 
such as Inlet line resistance (Rg) and Inertance (Lg) and discharge valve x*esistance 
(Rl), can also be easily Interpreted. 


The Inducer test loop of figure 16 can be simplified by assuming that: a 
large supply tank exists at location 8; a lax*ge discharge tank exists at location 2; 
all compliances, except that due to cavitation, are negligible; and the discharge 
tank Is sufficiently close to the Inducer discharge to make any Inertance between 
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tho two negligible. The resulting simplified model Is shown In flguro 90 and can 
be described by: 


1. Large supply tank maintained at constant prossuro (pg) 

2. Rigid Inlet line having resistance (Hg) and Inertance (Lg) 


3. Cavitation cavity at Inducer Inlet having a total volume (V c ) 
that Is dependent upon the Instantaneous values of Inlet static 
pressure (p 9 ) and Inlet mass flowrate (Wg) 

4, Inducer (or pump) turning at constant speed creating discharge 
static pressure (p^o) that Is dependent upon the Instantaneous 
values of Inlet static pressure (p 9 ), Inlet flowrate (Wo), and 
discharge flowrate (W 9 ) 

5, Rigid discharge line of negligible Inertance containing a dis- 
charge valve with resistance (Rj) 

6. Large discharge tank maintained at constant pressure (p 2 ). 


The linearized equations describing this system can be obtained from 
equations 12, 13, — With the simplifying assumptions these become: 

Inlet line: Ap 9 = - Rg AWg - Lg AWg 

Continuity: AW g - AW 9 = - p.AV c 

Cavitation: - p AV C = C p Ap 9 + C w AWg 

Performance: Ap 10 = G p Ap g + G w AW g - R p AW g 

Discharge line: ApjQ^RjAWg 
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The eight constants In the above equations aro represented graphically In 
figure 91, All constants are obtained from the known geometry and operating 
point; the steady-state head-flow, head falloff, and cavity volume maps; and the 
relationships of equations (17 and 27). Figure 91(a) Is the conventional inducer 
hoad falloff map. Figure 91(b) Is constructed from figure 91(a) by first plotting 
head falloff ratio (X) vs stondy-state flowrate (W) for the operating point Inlet 
pressure (p 9 ) and multiplying by noncnvltatlng pressuro rise (P rnc ). Figure 91(c) 
Is constructed by multiplying the noncnvltatlng head-flow map (P rac ) by the opera- 
ting point head falloff ratio (\). Figures 91(d) and 91(e) are tho predicted cavity 
volume maps. Figures 91(f) and 91(g) are predicted or moasurod static pressure 
drop curves for the Inlet and discharge linos, respectively. Figure 91(b) Is con- 
structed from the known geometry of the Inlet lino. The terms (U t $/A.) and 
(U t 0/A l )(A l /A d )2 in figures 91(b) and 91(c) result from the velocity heads at 
Inducer Inlet and discharge and the fact that the maps represent Inducer total- 
to-total pressure rise while It Is static pressure at Inducer Inlet and discharge 
that is used In the equations describing the system. 

The five simplified system equations can be solved simultaneously to 
obtain the equation of motion for a typical parameter, such as cavity volume (V c ): 

|(R p + RjlLgCpf AV C + {(Rp + RxHRgCp - C w ) + G p Lg| AV C + 


{(Rp + R 1 )+R 8 G p -G w (AV c = 0 


The condition for a cavttatlon-induced instability is easily seen through 
analogy with the spring-mass-damper system of figure 92. The equation of 
motion for such a system is : 

m3? + B& + KX = 0 


The following analogy between the Inducer system and the spring-mass-damper 
system can be made: 


Mass : M~(Rp + R^LgCp 

Damper: B ~ (R p + RjMRgCp - C w ) + G p L 8 

Spring: K ~ (R p + R t ) + R 8 G p - G w 

Either system will be unstable If the damper or the spring constant becomes 
negative. The model then indicates that a cavitation-induced Instability will 
occur when: 


c w > R 8 C p + r£ + 8 R t 


or when: 


®w > R 8®p + R p + R 1 


(32) 

(33) 


Equation (32), Involving cavity volume changes, represents the "continuity” 
mechanism of an Instability, while equation (33), Involving Inducer performance, 
represents the "performance" mechanism mentioned' In paragraph 3.1. 
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Above head breakdown G p « 1 and G w approaches zero. Equation (33) 
does not predict an Instability, since (R p + Rj) and Rg are positive numbers. 
However equation (32) predicts that an Instability will occur If: 


: w > R 8 Cp + R . Ri 


(34) 


Both flow compliance (C w ) and pressure compliance (C p ) are zero under non- 
cavltatlng conditions. Since Inlet line inertance (Lg) Is a positive number, we 
see that noncavltatlng flow Is predicted to be stable, a necessary condition for 
an Instability model. As cavitation number Is lowered, both C w and C p Increase 
(always positively). An Instability starts when flow compliance (C w ) reaches a 
value that satisfies equation (34). When cavitation number Is lowered to the 
head breakdown point, pressure gain (G p ) takes on a large positive value. 
Referring to equation (32), this Is seen to be a stabilizing Influence and explains 
why oscillations have generally been noted to disappear when head breakdown 
begins. 


In the head breakdown region, flow gain (G w ) can also take on large posi- 
tive values. An Instability related to Inducer performance will start If G w 
reaches a value that satisfies equation (33). This type of Instability occurred In 
our tests, making It Impossible to set a steady-state operating point. 


Frequency of oscillation at the Initial Instant of the Instability can be 
estimated as the system undamped natural frequency. By analogy to the spring- 
mass-damper system the natural frequency Is: 


u> * s/ K/M = 


(Rp * R l> + R 8 G p - G w 
(Rp + R^LgCp 


Since the operating point of Interest Is usually above head breakdown 
(G p = 1 and G w approaches zero) and since Inlet line resistance Is usually 
negligible (Rg = 0), this reduces to: 

W = y/ 1/IgCp (35) 
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.u x „ Com P nr ‘ n K Ration (35) with oquntlon (34), with H a - 0, It can bo aeon 
that for a given pump (Imown H p ) and given hydraulic syatorn (known l.g and Hi), 
flow compliance (C’ w ) dote rm Inca the operating point at which the Inatablllty 
marts, and preaaure eompllaneo (C’ p ) determines the oscillation frequency, 
thus, It la neeoaaary to define Inductor cavitation In torms of both flowrate and 
pressure. 


'the simplified system model qualitatively explains the occurrence of 
cavitation-induced Instabilities and demonstrates that Inducer cavitation and the 
hydraulic systom In which the Inducer operates are linked together to create the 
conditions under which an Instability can occur. Both must be considered In an 
analysis of stability. The link also Implies that stability can bo achieved through 
changos to the hydraulic system as well as through changes to the Inducor. 
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SECTION (l 

roNCTJIHIONB AND RECOMMENDATIONS 

ri. j coNci.rmoNH 

(1) A system model etm explain tho occurrence of solf-lnduced Instiblllty In 
cnvitntlng Inducers and cun bn used to doflno stabilizing system changes, 

Tho model requires accurate definitions of Inducer and system parameters 
In terms of resistance, Inerlnnee, and compliance. Frequency prodlellonH 
are especially senHltlvo to pressure compliance. 

(2) Two mechanisms can cause Inducer Instabilities: (1) a "continuity" mech- 
anism that Is oporablo whenever flow-sensltlvo Inducer cavitation Is present; 
and (2) a "performance" mechanism that Is usually oporablo only In tho 
hoad breakdown region. Continuity Instabilities are of primary Interest 
becauso they can occur In tho usual Inducer operating region. Both in- 
stability mechanisms wore observed and correlated satisfactorily with 
predictions In this program. Since blade suction surface cavitation Is 
known to occur at Inlet pressures woll above those at which performance 

Is affected, and this cavitation Is Inherently sonsltlvo to flow perturbations, 
such cavitation Is believed to be primarily responsible for "continuity" 
type Instabilities. Tip clearance cavitation would contribute to "continuity" 
Instabilities but the magnitude of the contribution Is unknown. Tip clear- 
ance cavitation may also be stabilizing In overall effect If Its sensitivity 
to flow Is small relative to Its sensitivity to pressure. This can only be 
proved by developing a tip vortex cavity model. 

(3) Predictions of conditions under which an Instability may occur, the 
frequency of oscillations, and trends of the Instability with changing con- 
ditions can be predicted with a linearized analysis. Predictions of ampli- 
tude and frequency after the start of oscillations would require an analog 
treatment of the nonlinear relationships. Both linear and analog analyses 
should be conducted to thoroughly define the Instability, with the linear 
analysis being used to gain an appreciation of significant parameters and 
trends and the analog analysis to predict exact amplitudes and frequencies. 

(4) Cascade Instabilities are also explainable through the use of a system 
model, but tho test data generated In this program were Inconclusive 

In that regard. The "performance" mechanism Is believed to be the only 
operable mechanism In a two-dimensional cascade because Incidence 
angle, which Is analogous to Inducer flow coefficient, would be constrained 
to a constant value, rendering the "continuity" mechanism Inoporable. 
"Performance" type cascade Instabilities roqulro an Inflection In tho cas- 
cade loss vs Inlet pressure relationship, and such an Inflection was moas- 
urod In this program. 

0. 2 RECOMMENDA TIONS 

(1) A Unoarlzed model was used to generate predictions of operating con- 
ditions under which an Instability would be oxpoctod and the trends of 
tho Instability with changes In the system and operating conditions. Tho 
model predicted the onset of oscillations at high cavitation numbers and 
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Ililn was mibplantlalml by lent (lulu. Tho moimurnd it ni| >1 i ttirlon worn 
relatively small, howovor, and did not Ineroimn appreciably until ravh 
fntlon number wan lowered a miliHtanllal ammmt, Thene low amplitudes 
may bo quickly damped mil anil not ommo any pressure oscillations In an 
aotual engine system. It would bo desirable, therefore, to generate 
Instability frequency and amplltudo predictions over Iho lull range of 
cavitation number through Iho uho of an analog modol. Wo recommend 
that such an analog modol of Iho luduooi' system bo proparod and Hh 
predictions correlated with tool data to further substantiate Iho modeling 
technique. 


( 2 ) Hxpor l mental data that are available for substantiation of inducer blade 
suction surfaco cavity volume prodlctlons aro llmltod. Instability pro- 
dictions aro sonsltlvo to tho rolatlonshlp of cavity volume to flowrate. 
Additional Inducer tests and cavity modol refinement effort are rocom- 
mondod to provide a more accurate, substantiated definition of suction 
surfaco cavity volunio as a function of cavitation number and flow coef- 
ficient. 

(3) Tip clearance cavitation Influences Instabilities because of volume sensi- 
tivity to pressure and flowrate changes, but the magnitude of Its effect 

Is currently Indetermlnato. Analytical and experimental efforts to define 
thoso relationships Is recommended. 

(4) Detailed model predictions have only been generated for tho tost loop 
system, and Its general applicability to other systems and fluid should be 
demonstrated. At least one additional system In which Instabilities have 
been measured should be modeled and the predictions correlated with 
test results. 

(5) Several system changes that appear to have a stabilizing Influence have 
been Identified In this program. These changes are tangential Inlet flow 
Injection, prewhirl Inducing guide vanes, and added Inducer Intet pres- 
sure compliance. Such changes should be further defined analytically 
and demonstrated systematically In a ^basellnG*’ system whose character- 
istics are well known* The objectives would be (1) to demonstrate that 
the changes resulted In the predicted stability trends and (2) to provide 
preliminary definition of Inducer and/or system design changes that can 
be applied to obtain stability* 
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APPENDIX 

NOMENCLATURE 


Symbol 

A 

A 


Definition 

Cross-sectional area normal to cascade nozzle 
or inducer centerline 

Cascade nozzle moan cross sectional area 



c 

D 


f l» f 2* f 3 


G P 

G w 

H 

HI 

H 


t 

K 

k 

k 


1 


Damping coefficient 

Length of cavity collapse region 

Compliance of inducer loop water and pipes 

Cavitation pressure compliance 

Cavitation flow compliance 

Cascade noncavitating total pressure loss 
coefficient (C nc = 

Length of cavity to point of maximum height 

Cascade total pressure loss 

Functions representing inducer head falloff, head 
vs flow, and cavity volume maps 

Pressure gain 

Flow gain 

Total head 

Ideal Head 

Total -to-total head rise 

Maximum cavity height, step height, or 
static head 

Fluid Incidence angle 

Spring constant 

Local or cascade cavitation number 

Inducer cavitation number based on average 
static pressure and tip relative velocity 

Total length of cavity 


L 


Unit 

2 


L 


2 


FTL" 2 

L 

MF -1 L 2 

MF _1 L 2 

MM _1 T 

dimensionless 

L 


dimensionless 

fl~ 2 m‘ 1 t 

L 

L 

L 

L 

rad (deg) 

FL" 1 

dimensionless 

dimensionless 

L 
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Symbol 

Definition 

Unit 

A 

Cascade nozzle length 

L 

L 

Fluid Inertance of cascade inlet nozzle or 
Inducer loop 

FI,- 2 M-!t 2 

M 

Mass 

M 

N 

Rotative speed 

T“1 

p r 

Inducer total-to-total pressure rise 

FL“ 2 

P 

Static pressure 

fl “ 2 

Q 

Volume flowrate 

l^t " 1 

q 

Fluid velocity pressure 

fl ” 2 

R l. . . 10 

Static-to-static resistance of inducer loop 

fl" 2 m _1 t 

R 

C 

Cascade internal resistance 

fl _ 2 m _1 t 

b p 

Inducer internal resistance 

fl" 2 m“ 1 t 


Cascade nozzle resistance 

fl - 2 m" 1 t 

s 

Laplace variable 

— 

s 

Inducer or cascade tahgential blade spacing 

L 

U t 

Inducer tip speed 

lt " 1 

V 

Fluid absolute velocity (Assumed axial at 
inducer inlet) 

lt ” 1 

V 

c 

Total cavity volume 

L 3 

w 

Mass flowrate 

mt ” 1 

X 

Coordinate along cascade nozzle centerline or 
of oscillating mass 

L 

0 

Angle of characteristic equation root In the 
complex plane f or flow angle from axial 

rad (deg) 

a w 

Angle of cavity collapse 

rad (deg) 

A 

Denotes linear perturbation about steady-state 
operating point 



Inducer inlet flow coefDcient 

<* - V«t> 

dimensionless 


Damping ratio 

dimensionless 
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Symbol 

Definition 

Unit 

u> 

Undamped natural frequency 

T -1 

P 

Fluid density 

ml’ 3 

r 

Inducer normal blade spacing or step 
separation passage height 

L 

X 

Ratio of inducer head rise to noncavitating head 
rise or cascade pressure loss to noncavitating 
pressure loss 

dimensionless 


Inducer average head rise coefficient 
(t = gH r /U t 2 ) 

dimensionless 

Superscripts: 


— 

Average value 


t 

Relative to inducer 


Subscripts (unless otherwise indicated): 


d 

Condition at cascade or inducer discharge 

i v 

h 

Condition at inducer hub 


i 

Condition at cascade or inducer inlet 


n 

Condition at cascade nozzle inlet 


nc 

Noncavitating value 


t 

Condition above free surface in cascade plenum tank or 
condition at inducer tip 


u 

Component in tangential direction 


V 

Vapor pressure 


z 

Component in axial direction 


1,2,3... 

1.0 Condition at station in inducer loop 

\ 
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